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Constantly accumulating observational data continue to confirm that about 70 % of the energy 
density today consists of dark energy responsible for the accelerated expansion of the Universe. 
We present recent observational bounds on dark energy constrained by the type la supernovae, 
cosmic microwave background, and baryon acoustic oscillations. We review a number of theoretical 
approaches that have been adopted so far to explain the origin of dark energy. This includes the 
cosmological constant, modified matter models (such as quintessence, k-essence, coupled dark energy, 
unified models of dark energy and dark matter), modified gravity models (such as f{R) gravity, 
scalar-tensor theories, braneworlds) , and inhomogeneous models. We also discuss observational 
and experimental constraints on those models and clarify which models are favored or ruled out in 
current observations. 
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I. INTRODUCTION 

The discovery of the late-time cosmic acceleration reported in 1998 [l|, 01 based on the type la Supernovae (SN la) 
observations opened up a new field of research in cosmology. The source for this acceleration, dubbed dark energy 
Q, has been still a mystery in spite of tremendous efforts to understand its origin over the last decade 0-113 • Dark 
energy is distinguished from ordinary matter in that it has a negative pressure whose equation of state wn v, is close 
to —1. Independent observational data such as SN la (Tsl - fisj . Cosmic Microwave Background (CMB) [l9l - l22j - and 
Baryon acoustic oscillations (BAO) p3l - l25j | have continued to confirm that about 70 % of the energy density of the 
present Universe consists of dark energy. 

The simplest candidate for dark energy is the so-called cosmological constant A whose equation of state is wde = — 1- 
If the cosmological constant originates from a vacuum energy of particle physics, its energy scale is significantly larger 
than the dark energy density today [26] {p'^l ~ lO""^ GeV). Hence we need to find a mechanism to obtain the tiny 
value of A consistent with observations. A lot of efforts have been made in this direction under the framework of 
particle physics. For example, the recent development of string theory shows that it is possible to construct de Sitter 
vacua by compactifying extra dimensions in the presence of fluxes with an account of non-pcrturbative corrections 

The first step toward understanding the property of dark energy is to clarify whether it is a simple cosmological 
constant or it originates from other sources that dynamically change in time. The dynamical dark energy models 
can be distinguished from the cosmological constant by considering the evolution of wde- The scalar field models 
of DE such as quintessence [2^ - [38j and k-essence [39l - l4l| predict a wide variety of variations of wbe, but still the 
current observational data are not sufhcient to provide some preference of such models over the A-Cold-Dark-Matter 
(ACDM) model. Moreover, the field potentials need to be sufficiently flat such that the field evolves slowly enough to 
drive the present cosmic acceleration. This demands that the field mass is extremely small (m^ ~ lO^'^^eV) relative 
to typical mass scales appearing in particle physics 0, . However it is not entirely hopeless to construct viable 
scalar-field dark energy models in the framework of particle physics. We note that there is another class of modified 
matter models based on perfect fluids-so-called (generalized) the Chaplygin gas model [13, . If these models are 
responsible for explaining the origin of dark matter as well as dark energy, then they are severely constrained from 
the matter power spectrum in galaxy clustering [46| . 

There exists another class of dynamical dark energy models that modify General Relativity. The models that 
belong to this class are f{R) gravity (47l. l49l - [5^ (/ is a function of the Ricci scalar R), scalar-tensor theories (53l - [57| . 
and Dvali, Gabadadze and Porrati (DGP) braneworld model [58]. The attractive feature of these models is that 
the cosmic acceleration can be realized without recourse to a dark energy component. If we modify gravity from 
General Relativity, however, there are stringent constraints coming from local gravity tests as well as a number of 
observational constraints such as large-scale structure (LSS) and CMB. Hence the restriction on modified gravity 
models is in general very tight compared to modified matter models. We shall construct viable modified gravity 
models and discuss their observational and experimental signatures. 

In addition to the above mentioned models, there are attempts to explain the cosmic acceleration without dark 
energy. One example is the void model in which an apparent accelerated expansion is induced by a large spatial 
inhomogeneity |5ft-t6S] . Another example is the so-called backreaction model in which the backreaction of spatial 
inhomogeneities on the Friedmann-Lemaitre-Robertson- Walker (FLRW) background is responsible for the real accel- 
eration [64| - |66| . We shall discuss these models as well. 

This review is organized as follows. In Sec.|lT]we provide recent observational constraints on dark energy obtained 
by SN la, CMB, and BAO data. In Sec. IIIII we review theoretical attempts to explain the origin of the cosmological 
constant consistent with the low energy scale of dark energy. In Sec. IIVI we discuss modified gravity models of dark 
energy-including quintessence, k-essence, coupled dark energy, and unified models of dark energy and dark matter. In 
SeclVlwe review modified gravity models and provide a number of ways to distinguish those models observationally 
from the ACDM model. Sec.|Vl]is devoted to the discussion about the cosmic acceleration without using dark energy. 
We conclude in Sec. IVIII 

We use units such that c = h = 1, where c is the speed of light and h is reduced Planck's constant. The 
gravitational constant G is related to the Planck mass rripi = 1.2211 x 10^^ GeV via G = l/^ipi and the reduced 
Planck mass ilfpi = 2.4357 x 10^* GeV via = 8t:G — 1/A^pi, respectively. We write the Hubble constant today 
as Hq = 100 /ikmsec^^ Mpc^^, where h describes the uncertainty on the value Hq. We use the metric signature 
{-,+,+,+). 
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II. OBSERVATIONAL CONSTRAINTS ON DARK ENERGY 

The late-time cosmic acceleration is supported by a number of independent observations-such as (i) supernovae 
observations, (ii) Cosmic Microwave Background (CMB), and (iii) Baryon acoustic oscillations (BAO). In this section 
we discuss observational constraints on the property of dark energy. 



A. Supernovae la observations 

In 1998 Riess et al. [l| and Perlmutter et al. et al. 0] independently reported the late-time cosmic acceleration 
by observing distant supernovae of type la (SN la) . The line-element describing a 4-dimensional homogeneous and 
isotropic Universe, which is called the FLRW space-time, is given by [63| 



,da;^d2;'^ 



dr2 



1 -Xr2 



+ r2(d6l2 + sin2 



(1) 



where a{t) is the scale factor with cosmic time t, and K = +1, —1,0 correspond to closed, open and flat geometries, 
respectively. The redshift z is defined by z = ao/a — 1, where ao = 1 is the scale factor today. 

In order to discuss the cosmological evolution in the low- redshift regime {z < C(l)), let us consider non-relativistic 
matter with energy density pm and dark energy with energy density poE and pressure Pde, satisfying the continuity 
equations 



Pm + iHp„i = , 

Pde + 3i7(pDE + -Pde) = , 



(2) 
(3) 



which correspond to the conservation of the energy-momentum tensor T^^ for each component (V where 
V represents a covariant derivative). Note that a dot represents a derivative with respect to t. The cosmological 
dynamics is known by solving the Einstein equations 

G^, = SttGT^, , (4) 

where G^i/ is the Einstein tensor. For the metric ([T]) the (00) component of the Einstein equations gives [g^I 



tt2 SttG, 
H = ——{Pv 



- Pue) 



K 
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where H = a/a is the Hubble parameter. We define the density parameters 



SirGpn 



which satisfy the relation 51.,„ + f^DE + — 1 from Ec^. ([5]). Integrating Eqs. ([2]) and (jS]), we obtain 



SttGpde 



K 



K 



p^ = p'^^){l + zf 



(0) 

PDE = Pde exp 



3(1 -t- u;de) j~ 

az 

l + z 



(5) 



(6) 



(7) 



where "0" represents the values today and wde = Pde / Pde is the equation of state of dark energy. Plugging these 
relations into Eq. (O, it follows that 



H\z) = Hl 



^-^±^A~z\+nf{l + zf 



(8) 



The expansion rate H{z) can be known observationally by measuring the luminosity distance di^{z) of SN la. The 
luminosity distance is defined by dr = Lgj (AttJ-), where Ls is the absolute luminosity of a source and J- is an observed 
flux. It is a textbook exercise 

[i,[li,|63 to derive dL{z) for the FLRW metric ([T]): 



dL{z) 



l + z 



: sinh 



(0) 



dz 

m 



(9) 
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n^^' sinh( 



where E{z) = H{z)/Ho. The function fxix) = 1/^^"'°^ 
K = +1), fnix) = X (for K = 0), and /^(x) = sinhx (for K 
dL{z) + £d~z/H{~z), i.e. 



H{z) 



_d_ 

dl 



l + z 



Sl^-'x) can be understood as fK{x) — sinx (for 
1). For the fiat case [K = 0), Eq. © reduces to 



(10) 



Hence the measurement of the luminosity distance dL{z) of SN la allows us to find the expansion history of the 
Universe for z < C(l). 

The luminosity distance is expressed in terms of an apparent magnitude m and an absolute magnitude M of an 
object, as 



m 



M = 5 logi 



10 pc 



(11) 



The absolute magnitude M at the peak of brightness is the same for any SN la under the assumption of standard 
candles, which is around M ~ —19 [Hl^]- The luminosity distance dL{z) is known from Eq. pT|) by observing the 
apparent magnitude m. The redshift z of an object is known by measuring the wavelength Aq of light relative to its 
wavelength A in the rest frame, i.e. z = Ao/A — 1. The observations of many SN la provide the dependence of the 
luminosity distance dj^ in terms of z. 

Expanding the function © around z = 0, it follows that 



dUz) = — 



1 

Ho 
1 



z+il- 



0{z^ 



^+4(1 



3wuEn'^l + nf^ z^ + o{z'') 



(12) 



where a prime represents a derivative with respect to z. Note that, in the second line, we have used Eq. ([8]). In the 

presence of dark energy (wde < and ripg > 0) the luminosity distance gets larger than that in the flat Universe 

without dark energy. For smaller (negative) wde and for larger ripg this tendency becomes more significant. The 

open Universe without dark energy can also give rise to a larger value of c?i(z), but the density parameter O^-* is 

constrained to be close to from the WMAP data (more precisely, —0.0175 < fJ^-* < 0.0085 [2l|). Hence, in the low 
redshift regime (z < 1), the luminosity distance in the open Universe is hardly different from that in the fiat Universe 
without dark energy. 

As we see from Eq. (I12p . the observational data in the high redshift regime (z > 0.5) allows us to confirm the 
presence of dark energy. The SN la data released by Riess et al. [l| and Perlmutter et al. [2] in 1998 in the redshift 
regime 0.2 < z < 0.8 showed that the luminosity distances of observed SN la tend to be larger than those predicted 
in the fiat Universe without dark energy. Assuming a flat Universe with a dark energy equation of state wde = ^1 
(i.e. the cosmological constant), Perlmutter et al. Q found that the cosmological constant is present at the 99 % 
confidence level. According to their analysis the density parameter of non-relativistic matter today was constrained 
to be d^} = 0.28t o;gg (68 % confidence level) in the flat universe with the cosmological constant. 

Over the past decade, more SN la data have been collected b y a number of high-redshift surveys-such as SuperNova 
Legacy Survey (SNLS) [ISjlHubble Space Telescope (HST) [T^.ll5i. and "Equation of State: SupErNovae trace Cosmic 
Expansion" (ESSENCE) [3, [13] survey. These data also confirmed that the Universe entered the epoch of cosmic 
acceleration after the matter-dominated epoch. If we allow the case in which dark energy is different from the 
cosmological constant (i.e. wde —1), then observational constraints on wde and Slj^g (or VL^) are not so stringent. 
In Fig. [T]we show the observational contours on (wDE,^^m^) for constant wde obtained from the "UnionOS" SN la 
data by Kowalski et al. [18]. Clearly the SN la data alone are not yet sufficient to place tight bounds on wde- 

In the flat Universe dominated by dark energy with constant wde, it follows from Eq. ([S]) that ~ iJp 57^^(1 + 
2,")3(i+i«de) Q-3(i+ii)DE)^ Integrating this equation, we flnd that the scale factor evolves as a c>c i2/(3(i+M>DE)) fQj. 
wde > —1 and a cx e^* for wde — ~1- The cosmic acceleration occurs for —1 < wde < —1/3. In f act, Fig. [T]shows 
that wde is constrained to be smaller than —1/3. If wde < — 1, which is called phantoms or ghosts |68l], the solution 
corresponding to the expanding Universe is given by a oc (ts — ^)2/(3(i+"'de))^ where tg is a constant. In this case the 
Universe ends ait — tg with a so-called big rip singularity [69|, |70j at which the curvature grows toward infinity.^ The 



^ There are other classes of finite-time singularities studied in Refs. [7ll479ll . In some cases quantum effects can moderate such singularities 
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UnionOS SNe 
UnionOS SNe w/ sys 




Figure 1: 68.3 %, 95.4 %, and 99.7 % confidence level contours on uide and (denoted as w and ilm in the figure) constrained 
by the UnionOS SN la datasets. The equation of state iude is assumed to be constant. From Ref. [l^ . 



current observations allow the possibility of the phantom equation of state. We note, however, that the dark energy 
equation of state smaller than —1 does not necessarily imply the appearance of the big rip singularity. In fact, in 
some of modifi ed g ravity models such as f{R) gravity, it is possible to realize wue < — 1 without having a future big 
rip singularity [83]. 

If the dark equation of state is not constant, we need to parametrize wde as a function of the redshift z. This 
smoothing process is required because the actual observational data have discrete values of redshifts with systematic 
and statistical errors. There are several ways of parametrizations proposed so far. In general one can write such 
parametrizations in the form 

■Wde{z) = ^ WnXn{z) , (13) 
ri=0 

where n's are integers. We show a number of examples for the expansions: 

(i) Redshift: a;„(z) = z", (14) 



(ii) Scale factor : x„(z) = (1 — a)" = (^Y~i~y ' ^^^"^ 

(iii) Logarithmic : Xn(z) — [ln(l + z)]" . (16) 
The paramctrization (i) was introduced by Huterer and Turner [l] and Weller and Albrecht ^ with rt < 1, i.e. 



wde =^ wq + wiz. Chevalier and Polarski and Lindcr [8u\ proposed the parametrization (ii) with n < 1, i.e. 

z 

wbe{z) ^ wq + wi{l - a) = wq + wi— — . (17) 

1 + z 

This has a dependence wde(^) = + wi for z — >■ oo and wbe{z) — )■ wq for 2; — )■ 0. A more general form, 
we)e{z) — wq + Wiz/{1 + zY, was proposed by Jassal et al. [s^]. The parametrization (iii) with n < 1 was introduced 
by Efstathiou 'ss]. A functional form that can be used for a fast transition of wde(-z) was also proposed [s^ [90| . 
In addition to the parametrization of wde, a number of authors assumed parametric forms of dL{z) [9l|, or H{z) 



92 



9411. Many works placed observational constraints on the property of dark energy by using such parametrizations 



[zllMH. 
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■ WMAP+Ho+SN 

■ WMAP+BAO:+Ho+SN 
Q WMAP+BAO:+H(,+Da,+SN 



Wo 

Figure 2: Observational constraints on the parameters wo and wi (denoted as Wa in the figure) for the parametrization (I17p . 
The contours show the 68 % and 95 % confidence level from WMAP+//o+SN (red), WMAP+BAO+i^o+SN (blue), and 
WMAP+BAO+Jfo+DAt+SN (black), for a flat universe. From Ref. [H. 



In Fig. [2] we show the SN la constraints combined with other measurements such as the WMAP 7-year 22| and 
the BAO data psj . The parametrization pT|) is used in this analysis. The Gaussian prior on the present-day Hubble 
constant 'llO], Hq — 74.2 ± 3.6 km sec~^ Mpc~^ (68 % confidence level), is also included in the analysis (obtained 
from the magnitude-redshift relation of 240 low-z SN l a at z < 0.1). In Fig. [21 "-Da*" means a constraint coming from 
the measurement of gravitational lensing time delays |lllj . The joint constraint from WMAP+BAO+ifo+-DAt+SN 
gives the bound 



Wo = -0.93 ±0.13, 



wi = -0.4lt°;^? , 



(18) 



at the 68 % confidence level. Hence the current observational data are consistent with the flat Universe in the presence 
of the cosmological constant {wq ~ —1,^1 — 0). 



B. CMB 



The temperature anisotropics in CMB are affected by the presence of dark energy. The position of the acoustic 
peaks in CMB anisotropics depends on the expansion history from the decoupling epoch to the present. Hence the 
presence of dark energy leads to the shift for the positions of acoustic peaks. There is also another effect called the 
Integrated-Sachs- Wolfe (ISW) effect |ll2j | induced by the variation of the gravitational potential during the epoch 
of the cosmic acceleration. Since the ISW effect is limited to large-scale perturbations, the former effect is typically 
more important. 

The cosmic inflation in the early Universe 113l4ll6j predicts nearly scale-invariant spectra of density perturbations 
through th e qua ntum fluctuation of a scalar fleld. This is consistent with the CMB temperature anisotropics observed 
by COBE [ll3 and WMAP The perturbations are "frozen" aft er the sca le A — {2tt /k)a {k is a comoving 

wavenumber) leaves the Hubble radius H^^ during inflation (A > H~^) [ll8l . Ill9| . After inflation, the perturbations 
cross inside the Hubble radius again (A < H^^) and they start to oscillate as sound waves. This second horizon 
crossing occurs earlier for larger k (i.e. for smaller scale perturbations). 

We define the sound horizon as rs{ri) = dfjCsifj), where Cs is the sound speed and di] = a^^dt. The sound speed 
squared is given by 



(19) 



where pt and pj are the energy densities of baryons and photons, respectively. The characteristic angle for the location 
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of CMB acoustic peaks is [l20j | 



(^dcc) 



(20) 



(c) 

where d)^ is the comoving angular diameter distance related with the luminosity distance via the duality relation 

d^-* = dL/{l + z) [l^, [63], and Zdoc — 1090 is the redshift at the decoupling epoch. The CMB multipole Ia that 
corresponds to the angle ||^D|) is 



s (^dcc) 



(21) 



Using Eq. ^ and the background equation 3i/^ = 8TrG{pm + Pr) for the redshi ft z > Zde c (where pm and pr are the 
energy density of non-relativistic matter and radiation, respectively) , we obtain jl2ll Il22l | 



In 



^/ Rs (adec) + Rs (acq) + ^1 + -Rs (adec) 
1+ ^i?,(aeq) 



7^, 



(22) 



where = ri[°^/i^ and ui^ = fl^^h'^, and TZ is the so-called CMB shift parameter defined by [12 



n 



(0) 



sinh 



dz 
E(z) 



(23) 



The quantity Rg = 3/9f,/(4/9-y) can be expressed as 



Rs(a) = (3ajf,/4a;^)a . 



(24) 



In Eq. (|22]) , adoc and Ooq correspond to the scale factor at the decoupling epoch and at the radiation-matter equality, 
respectively. 

The change of cosmic expansion history from the decoupling epoch to the present affects the CMB shift parameter, 
which gives rise to the shift for the multipole The general relation for all peaks and troughs of observed CMB 
anisotropies is given by |l24| 



(25) 



where m represents peak numbers (m = 1 for the first peak, m — 1.5 for the first trough,...) and <j)m is the shift of 
multipoles. For a given cosmic curvature JT^^the quantity (j)rn depends weakly on cjf, and w,„ = d^h"^. The shift of 
the first peak can be fitted as (jii = 0.265 Il24| . The WMAP 5-year bound on the CMB shift parameter is given by 

m 



n = 1.710 ±0.019, 



(26) 



at the 68 % confidence level. Taking TZ — 1.710 together with other values cot = 0.02265, uJm = 0.1369, and 
Lo-y = 2.469 X 10"^ constrained by the WMAP 5-year data, we obtain £a ^ 300 from Eq. ^2^. Using the relation (^5)) 



with (pi — 0.265 we find that the first acoustic peak corresponds to £i 

In the flat Universe {K — 0) the CMB shift parameter is simply given by TZ 



rim' (i.e. for larger 51de)' ^ tends to be smaller. For the cosmological constant (iude = ^1) the normalized Hubble 



l(0) 



220, as observed in CMB anisotropies. 

(0) rz, 
JQ 



' dz/E{z). For smaller 



expansion rate is given by E{z) = [fim'' (1 + z) 



-^de]^^^- Under the bound 



the density parameter is constrained 



to be 0.72 < Slpg < 0.77. This is consistent with the bound coming from the SN la data. One can also show that. 



for increasing wde, the observationally allowed values of ilm^ gets larger. However, TZ depends weakly on the wde- 
Hence the CMB data alone do not provide a tight constraint on wde- In Fig. [3] we show the joint observational 
constraints on wbe and Q,m^ (for constant wde) obtained from the WMAP 5-year data and the UnionOS SN la data 
[l8l |. The joint observational constraints provide much tighter bounds compared to the individual constraint from 
CMB and SN la. For the flat Universe Kowalski et al. 18] obtained the bounds wde = 



= 0.265 



n nc;c+0. 060-1-0. 059 i 
u.yoo_Q Qgg_Q QgQ anu 



+0.022+0.018 



(with statistical and systematic errors) from the combined data analysis of CMB and SN la. 



See also Refs. 9^ 1251 ^ 133 [ for related observational constraints 
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Figure 3: 68.3 %, 95.4 % and 99.7 % confidence level contours on woe and (denoted as w and Q,m in the figure, respectively) 
for a fiat Universe. The left panel illustrates the individual constraints from SN la, CMB, and BAO, as well as the combined 
constraints (filled gray contours, statistical errors only). The upper right panel shows the effect of including systematic errors. 
The lower right panel illustrates the impact of the Supernova Cosmology Project (SCP) Nearby 1999 data. From Ref. [l^ . 



C. BAO 



The detection of baryon acoustic oscillations first reported in 2005 by Eisenstein et. al. in a spectroscopic 
sample of 46,748 luminous red galaxies observed by the Sloan Digital Sky Survey (SDSS) has provided another test 
for probing the property of dark energy. Since baryons are strongly coupled to photons prior to the decoupling epoch, 
the oscillation of sound waves is imprinted in baryon perturbations as well as CMB anisotropics. 

The sound horizon at which baryons were released from the Compton drag of photons determines the location of 
baryon acoustic oscillations. This epoch, called the drag epoch, occurs at the redshift Zd- The sound horizon at z = 
is given by r sizd) = J^'' dry 0^(77), where Cs is the sound speed. According to the fitting formula of Zd by Eisenstein 
and Hu |134| . Zd and rs{zd) are constrained to be around Zd ~ 1020 and rs{zd) ~ 150 Mpc. 

We observe the angular and redshift distributions of galaxies as a power spectrum P{k±, fc||) in the redshift space, 
where k± and are the wavenumbers perpe ndicular and parallel to the direction of light respectively. In principle 
we can measure the following two ratios [135l | 

a ( \ ^'^(^'i) X ( \ r,{zd)H{z) 

dy{z) c 

(c) . 

where the speed of light c is recovered for clarity. In the first equation d\ is the comoving angular diameter distance 

related with the proper angular diameter distance d-A via the relation d^^ = dA/a = ^^(1 + z). The quantity 9s{z) 
characterizes the angle orthogonal to the line of sight, whereas the quantity 5zs corresponds to the oscillations along 
the line of sight. 

The current BAO observations are not sufficient to measure both 9s{z) and dzJz) independently. From the 
spherically averaged spectrum one can find a combined distance scale ratio given by [135l | 

or, alternatively, the effective distance ratio [1^ 

Dviz) = [(1 + z)24(z)cz/iJ(z)]'/' . (29) 

In 2005 Eisenstein et al. Q obtained the constraint Dv{z) = 1370 ± 64 Mpc at the redshift z = 0.35. In 2007 
Percival et al. p^l measured the effective distance ratio defined by 



rBAo{z) = rs{zd)/Dv{z) , 



(30) 
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at the two redshifts: rBAo{z = 0.2) = 0.1980 ± 0.0058 and rBAo(^ = 0.35) = 0.1094 ± 0.0033. This is based on the 
data from the 2-degree Field (2dF) Galaxy Redshift Survey. These data provide the observational contour of BAO 
plotted in Fig. [3l From the joint data analysis of SN la WMAP 5-year [lH, and BAO data [H, Kowalski et 

al. [13 placed the constraints wde = -0.969i;];{]^|^(stat)+5];[5^^(sys) and VL^^ = 0.2741J] ;]}^(stat)tg;[]i^(sys) for the 
constant equation of state of dark energy. 

The recent measurement of the 2dF as well as the SDSS data provided the effective distance ratio to be tbao {z — 
0.2) = 0.1905±0.0061 and rBA_o{z = 0.35) = 0.1097 ± 0.0036 ^2?!. Using these data together with the WMAP 7-year 
data [22] and the Gaussian prior on the Hubble constant Hq = 74.2 ± 3.6 km sec~^ Mpc~^ jllQ] . Komatsu et al. 
[2^ 1 derived the constraint wde = —1.10 ± 0.14 (68 % confidence level) for the constant equation of state in the flat 
Universe. Adding the high-z SN la in their analysis they found the most stringent bound: wde = — 0.980 ± 0.053 (68 
% confidence level). Hence the ACDM model is well consistent with a number of independent observational data. 



Finally we should mention t hat there are other constrain ts coming from the cosmic age [136| , large-scale clustering 
Il37l - ll3a . gamma ray bursts [T40l - [l4^ . and weak lensing [l45l - [l50ll . So far we have not found strong evidence for 



supporting dynamical dark energy models over the ACDM model, but future high-precision observations may break 
this degeneracy. 

III. COSMOLOGICAL CONSTANT 

The cosmological constant A is one of the simplest candidates of dark energy, and as we have seen in the previous 
section, it is favored by a number of observations. However, if the origin of the cosmological constant is a vacuum 
energy, it suffers from a serious problem of its energy scale relative to the dark energy density today ^265. The 
zero-point energy of some field of mass m with momentum k and frequency uj is given hy E = uj /2 = \/fc^^T?/2. 
Summing over the zero-point energies of this field up to a cut-off scale fcmax (^ fn)^ we obtain the vacuum energy 
density 



(27r)3 2 

Since the integral is dominated by the mode with large k (3> m), we find that 

" 47rfc2dfc 1 . k: 



y/k"^ + m2 . (31) 



(2^)3 2^" 16^2 • ^ 

Taking the cut-off scale fcmax to be the Planck mass mpi, the vacuum energy density can be estimated as pvac — 

10''* GeV*. This is about 10"^^-^ times larger than the observed value p^g ~ 10"'''' GeV''. 

Before the observational discovery of dark energy in 1998, most people believed that the cosmological constant is 
exactly zero and tried to explain why it is so. The vanishing of a constant may imply the existence of some symmetry. 
In supersymmetric theories the bosonic degree of freedom has its Fermi counter part which contributes to the zero 
point energy with an opposite sign^. If supersymmetry is unbroken, an equal number of bosonic and fermionic degrees 
of freedom is present such that the total vacuum energy vanishes. However it is known that supersymmetry is broken 
at sufficient high energies (for the typical scale Msusy ~ 10'^ GeV). Therefore the vacuum energy is generally non-zero 
in the world of broken supersymmetry. 

Even if supersymmetry is broken there is a hope to obtain a vanishing A or a tiny amount of A. In supergravity 
the ory the effective cosmological constant is given by an expectation value of the potential V for chiral scalar fields 



^p' [1511: 



y((^,V3*) = e^^ ^ D,W{K'^ ){DjW)* ~?,k'\WY , (33) 

where K and W are the so-called Kahler potential and the superpotential, respectively, which are the functions of 
and its complex conjugate ip" . The quantity K^^ is an inverse of the derivative Kij- = d'^Kjdip^dip^ , whereas the 
derivative D.^W is defined by D,W = dW/dtp^ + K'^W{dK/dtp'). 

The condition DiW 7^ corresponds to the breaking of supersymmetry. In this case it is possible to find scalar 
field values leading to the vanishing potential {V — 0), but this is not in general an equilibrium point of the potential 



^ The readers who are not familiar with supersymmetric theories may consult the books Il5ll . Il52(l . 
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V. Nevertheless there is a class of Kahler potentials and superpotentials giving a stationar y sca lar-field configuration 
atV = 0. The gluino condensation model in Eg, x Eg superstring theory proposed by Dine jl53j | belongs to this class. 
The reduction of the 10-dimensional action to the 4-dimensional action gives rise to a so-called modulus field T. This 
field characterizes the scale of the compactified 6-dimensional manifold. Generally one has another complex scalar 
field S corresponding to 4-dimensional dilaton/axion fields. The fields T and S are governed by the Kahler potential 

K{T, S) = -(3/k2) In {T + T*) - {l/K^)\xi [S + S*) , (34) 

where [T + T*) and {S + S*) are positive definite. The field S couples to the gauge fields, while T does not. An 
effec tive superpotential for S can be obtained by integrating out the gauge fields under the use of the i?-invariance 
[l54j : 

W{S) = Ml, [ci + C2 exp(-35/2c3)] , (35) 

where ci,C2, and C3 are constants. 

Substituting Eqs. ([M)) and ([55)1 into Eqs. ([55)1 . we obtain the field potential 

y - (r-^T-)3(^ + ^.) (^^^)^^"(^^^)* 



^4 



{T + T*Y{S + S*) 



ci + C2exp(-35/2c3) <! 1 + -^(5 + 5* 

2c3 



(36) 



where, in the first line, we have used the property {DtW)K^^ (DtW)* — 3k'^\W\'^ for the modulus term. This 
potential is positive because of the cancellation of the last term in Eq. (|33|) . The stationary field configuration with 
V = is reahzed under the condition DgW = dW/dS - W/{S + S*) = 0. The derivative, DtW = K^WdK/dT = 
—3W/ (T-l-T*), does not necessarily vanish. When DtW 7^ the supersymmetry is broken with a vanishing potential 
energy. Therefore it is possible to obtain a stationary field configuration with ^ = even if supersymmetry is broken. 

The discussion above is based on the lowest-order perturbation theory. This picture is not necessarily valid to 
all finite orders of perturbation theory because the non-supersymmetric field configuration is not protected by any 
symmetry. Moreover some non-perturbative effect can provide a large contribution to the effective cosmological 
constant [43!]. The so-called flux compactification in type IIB string theory allows us to realize a metastable de 
Sitter (dS) vacuum by taking into account a non-perturbative correction to the superpotential (coming from brane 
instantons) as well as a number of anti D3-branes in a warped geometry Hence it is not hopeless to obtain a 
small value of A or a vanishing A even in the presence of some non-perturbative corrections. 

Kachru, Kallosh, Linde and Trivedi (KKLT) [23| constructed dS solutions in type II string theory compactified on 
a Calabi-Yau manifold in the presence of flux. The construction of the dS vacua in the KKLT scenario consists of 
two steps. The flrst step is to freeze all moduli fields in the fiux compactification at a supersymmetric Anti de Sitter 
(AdS) vacuum. Then a small number of the anti D3-brane is added in a warped geometry with a throat, so that 
the AdS minimum is uplifted to yield a dS vacuum with broken supersymmetry. If we want to use the KKLT dS 
minimum derived above for the present cosmic acceleration, we require that the potential energy Vds at the minimum 
is of the order of Vds — 10"'''^ GeV"^ . De pending on the number of fluxes there are a vast of dS vacua, which opened 
up a notion called string landscape (l55j |. 



The question why the vacuum we live in has a very sm all energy density among many possible vacua has been 



a very sm 

sometimes answered with the a nthr opic principle [I56lll57j . Using the anthropic arguments, Weinberg put the bound 



on the vacuum energy density [158 1 

- 10-123 m^j < PA < 3 X lO-i^^m^i . (37) 

The upper bound comes from the requirement that the vacuum energy does not dominate over the matter density 
for the redshift z > 1. Meanwhile the lower bound comes from the condition that pA does not cancel the present 
cosmological density. Some people have studied landscape statistics by considering the relative abundance of long- 
lived low-energy vacua satisfying the bound p7p |l59l4l62| . These statistical approaches are still under study, but it 
will be interesting to pursue the possibility to obtain high probabilities for the appearance of low-energy vacua. 

Even in 1980s there were some pioneering works for finding a mechanism to make the effective cosmological constant 
small. For example, let us consider a 4- form field F^^'^^'^ expressed by a unit totally anti-symmetric tensor e'^'^^°' , as 
p,,uXa ^ pgMi^A^T constant). Then the energy density of the 4-form field is given by Ff""^" F^^^xa / {'2. ■ 4!) = 0^/2. 

Taking into account a scalar field </> with a potential energy F(0), the total energy density is A = V{(t)) + c? /2. In 
1984 Linde |163{ considered the quantum creation of the Universe and claimed that the final value of A can appear 
with approximately the same probability because V{<j)) can take any initial value such that A k, m^,. 
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In 1987-1988 Brown and Teltelboim jl64l . Il65l | studied the quantum creation of closed membranes by totally an- 
tisymmetric tensor and gravitational fields to neutralize the effective cosmological constant with small values. The 
constant c appearing in the energy density of the 4-form field can be quantized in integer multiples of the membrane 
charge g, i.e. c = nq. If we consider a negative bare cosmological constant — Af, (as in the KKLT model) in the presence 
of the flux energy density then the effective gravitational constant is given by A = — Ab -I- The field 

strength of the 4-form field is slowly discharged by a quantum Schwinger pair creation of field sources [nq — > (n — l)g]. 
However, in order to get a tiny value of A consistent with the dark energy density today, the membrane change q is 
constrained to be very small (for natu ral choices of Ab) [166]. 

In 2000 Bousso and Polchinski jl67l | considered multiple 4-form fields that arise in M-theory compactifications and 
showed that the small value of A can be explained for natural choices of q. More precisely, if we consider J(> 1) 
4-form fields as well as J membrane species with charges qi, q2, - ■ ■ ,qj and the quantized flux p^'-'^^'^ = riiqit^'^^'^ , the 
effective cosmological constant is given by 

J 

A = -A, + Y,nh'/2- (38) 

1=1 

Bousso and Polchinski 167] showed that, for natural values of charges {qi < 0(0.1)) there exists integers Ui such that 
2Af, < J2Li ^hf < 2(A6 + AA) with AA « 10-^7 GeV*. This can be realized for J > 100 and Af, « m^j. 

There are some interesting works for decoupling A from gravity. In the cascading gravity scenario proposed in 
Ref. [168i] the cos molo gical constant can be made gravitationally inactive by shutting off large-scale density pertur- 
bations. InRef. |l69l | an incompressible gravitational Aether fluid was introduced to degraviate the vacuum. In 

Refs. [Mlinl 

Padmanabhan showed an example to gauge away the cosmological constant from gravity according to 
the variational principle different from the standard method. See also Refs. 172 184] for other possibilities to solve 
the cosmological constant problem. If the cosmological constant is completely decoupled from gravity, it is required 
to find alternative models of dark energy consistent with observations. 

In the subsequent sections we shall consider alternative models of dark energy, under the assumption that the 
cosmological constant problem is solved in such a way that it vanishes completely. 



IV. MODIFIED MATTER MODELS 



In this section we discuss "modified matter models" in which the energy-momentum tensor T^^, on the r.h.s. of the 
Einstein equations contains an exotic matter source with a negative pressure. The models that belong to this class 
are quintessence, k-essence, coupled dark energy, and generalized Chaplygin gas. 



Quintessence 



A canonical scalar field 4> responsible for dark energy is dubbed quintessence [3^ [37 
earlier works). The action of quintessence is described by 



S = 



d^a 



-9 



s 



M ■ 



(see also Refs. [28H3J for 



(39) 



where i? is a Ricci scalar, and is a scalar field with a potential V{4>). As a matter action Sm, we consider perfect 
fluids of radiation (energy density pr, equation of state Wr = 1/3) and non-rclativistic matter (energy density pm, 
equation of state — 0). 

In the flat FLRW background radiation and non-rclativistic matter satisfy the continuity equations pr + AH pr = 
and pm + iHpm = 0, respectively. The energy density p^ and the pressure of the fleld are p^ = 0^/2 + V{(j)) and 
= 0^/2 — V^(0), respectively. The continuity equation, p^ + 'iH{p^ + P^) — 0, translates to 



+ 3i/0 + 

where F ^ = AV/ d<f). The fleld equation of state is given by 



0, 



_ _ 0^ - 2^(0) 



(40) 



(41) 



12 



From the Einstein equations (jj]) we obtain the following equations 



H 



K 

y 



K 

Y 



Vict)) + Prr. 
4 



Pr 



(42) 
(43) 



Although {pr,Pm\ ^ P4> during radiation and matter eras, the field energy density needs to dominate at late 
times to be responsible for dark energy. The condition to realize the late-time cosmic acceleration corresponds to 
UJ0 < —1/3, i.e. (fP' < V{4)) from Eq. (HJ). This means that the scalar potential needs to be flat enough for the field 
to evolve slowly. If the dominant contribution to the energy density of the Universe is the slowly rolling scalar field 
satisfying the condition 0^ <C V{(f)), we obtain the approximate relations iHcf) + ~ and 2>H^ ~ n^V{4)) from 
Eqs. (j40| and (j42|). respectively. Hence the field equation of state in Eq. (|4ip is approximately given by 



-1 + 2e^/3, 



(44) 



where = {V,^/V) / (2k^) is the so-called slow-roll parameter |118| . During the accelerated expansion of the Universe, 
ts is much smaller than 1 because the potential is sufficiently flat. Unlike the cosmological constant, the field equation 
of state deviates from —1 [w^ > — !)■ 

Introducing the dimensionless variables Xi = ncj)/ {^/6H) X2 = k VV/{V3H), and X3 = Ky/p^/{V3H), we obtain 
the following equations from Eqs. (gO]), (gH), and (gS]) fol. Isa. [iSSi flSej : 



x[ — Sxi + -^y^x\ + 2^l("^ + "^^1 ^ "^^2 + ^3) : 
X2 = —-^^XiX2 + 2^2(3 + 'ix\ — ix\ + 0:3) , 



-2x3, + ^^3(3 + ix\ 



where a prime represents a derivative with respect to = In a, and A is defined by A 
parameters of the field, radiation, and non-relativistic matter are given by il^ = x\ - 



(45) 
(46) 
(47) 



= -V^4,I{kV). The density 



1 



2:3, respectively. One has constant A for the exponential potential [33| 



(48) 



in which case the fixed points of the system p5)) - (l47)) can be derived by setting x[ — Q {i — \, 2, 3). The fixed point 
that can be used for dark energy is given by 



(xi,X2,X3) = (a/V6, - A2/6,o) , u;0 = -1 + AV3, 



(49) 



The cosmic acceleration can be realized for < —1/3, i.e. A'^ < 2. One can show that in this case the accelerated 
fixed point is a stable attractor |35i]. Hence the solutions finally approach the fixed point after the matter era 
[characterized by the fixed point (xi, a;2, xa) = (0,0,0)]. 
If A varies with time, we have the following relation 



A' = -^/&\^{T -\)xi 



(50) 



where F = VV^^^/V"^^. For monotonically decreasing potentials one has A > and xi > for < and A < and 
xi <0 for V A, > 0. If the condition 



VI 



> 1. 



(51) 



is satisfied, the absolute value of A decreases toward irrespective of the signs of V^ [331 • Then the solutions finally 
appr oach the accelerated "instantaneous" fixed point (|49p even if A^ is larger than 2 during radiation and matter eras 
[l85l . [186]. In this case th e fie ld equation of state gradually decreases to —1, so the models showing this behavior are 
called "freezing" models flST]. The condition ^ is the so-called tracking condition under which the field density 
eventually catches up that of the background fluid. 
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Figure 4: The allowed region in the (ui^jUi^) plane for thawing and freezing models of quintessence (^10 is denoted as w in the 
figure). The thawing models correspond to the region between two curves: (a) w'^ — 3(1 + W0) and (b) w'^ = l + w^, whereas the 
freezing models are characterized by the region between two curves: (c) w'^ = 0.21^^(1 + w^) and (d) w'^ = iw^{l + w^). The 
dotted line show s the border between the acceleration and deceleration of the field (0 = 0), which corresponds to w'^ = 2>{1+'W^Y' . 
From Ref. [iSTII. 



A representative potential of the freezing model is the inverse power-law potential V{(t)) = Af^+" ^~" (n > 0) 
[31I [38| , which can appear in the fermion condensate model as a dynamical supersymmetry breaking |l88l | . In this 
case one has r = (n + l)/n>l and hence the tracking condition is satisfied. Unlike the cosmological constant, even if 
the field energy density is not negligible relative to the background fluid density around the beginning of the radiation 
era, the fleld eventually enters the tracking regime to lead to the late-time cosmic acceleration [SSs] • Another example 
of freezing models is V{(t)) = M^'^"(j)~" exp{a(fP / m'^i) , which has a minimum with a positive energy density at which 
the field is eventually trapped. This potential is motivated in the framework of su pergravity [l89l |. 

There is another class of quintessence potentials called "thawing" models jl87l |. In thawing models the field with 
mass has been frozen by the Hubble friction (i.e. the term Hcj)) until recently and then it begins to evolve 
after H drops below m^. At early times the equation of state of dark energy is Wcf, ~ —1, but it begins to grow 
for H < m^. The representative potentials that belong to this class are (a) V{(j)) = + M ^~"^ " (n > 0) and 
(b) V{(j)) = cos^ {(j)/ f) . The potential (a) with n = 1 wa s or iginally proposed by Linde |l9C| to replace the 
cosmological constant by a slowly evolving scalar field. In Ref. |l91i | this was revised to allow for negative values of 
V{(p). The universe will collapse in the future if the system enters the region with V{(j)) < 0. T he potential (b) is 
motivated by the Pseudo-Nambu-Goldstone Boson (PNGB), which was introduced in Ref. [192] in response to the 
first tentative suggestions for the existence of the cosmological constant. The small mass of the PNGB model required 
for dark energy is protected against radiative corrections, so this model is favored theoretically. In fact there are a 
number of interesting wo rks to ex plain the small energy scale M w 10~^ eV required for the PNGB quintessence in 
supersymmetric theories |l93l - ll96| . See Refs. [197 202, 204 208] for the construction of quintessence potentials in the 
framework of supersymmetric theories. 

The freezing models and the thawing models are characterized by the conditions w'^ = dw^/dN < and 
w'^ > 0, respectively. More precisely the allowed regions for the freezing and thawing models are given by 
3w^{l + w^) <w'^< 0.2w^{l + w^) and 1 + ^ ^'4, ^ 3(1-1- w;^), respectively |l87l | (see Ref. [209] for details). These 
regions are illustrated in Fig. U) While the observational data up to now are not sufficient to distinguish freezing and 
thawing models by the variation of w^, we may be able to do so with the next decade high-precision observations. 

There is also another useful measure calle d "statefinder" j210| by which the quintessence models can be distinguished 
from the ACDM model (see also Ref. [21l| ). The statefinder parameters are defined by 

- aH^ di3 ' 3{q - 1/2) ' ' 



where q — 



~a/{aH'^) is the deceleration parameter. Let us consider the case in which the dark energy density pde 
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satisfies the continuity equation pde + 3i?(l + u'de)pde = 0. In Einstein gravity in which the Friedmann equation 
3H^ = K^ipuE + Pm) holds, we obtain 

9wDEf^DE , , WdE 

r = l^ s, s = 1 + wde-— , (53) 

2 3HwuE 

where f^DE = h'^Pde/{3H'^). The ACDM model (wde — —1) corresponds to a point (r, s) = (1,0) in the (r, s) plane, 

but the quintessence models are characterized by the curves in the region r < 1 and s > 93]. The Chaplygin gas 

model we discuss in Sec. IIVDI gives rise to the curves in the region r > 1 and s < 0. Hence one can distinguish 

between dark energy models by using the statefinders defined in Eq. (|52p . 

It is possible to reconstruct quintessence potentials from the observational data of SN la. Neglecting the contribution 

of radiation and using the relat ion P m = Pm\^ + z)^ for non-relativistic matter, we get the following equations from 

Eqs. dm and dlSl) [9ll[9l [2T2l . [2T3| : 



fd(j)Y 1 din E{z) 3fi^^ 1 + z 



2 \dz J 1 + z dz 2 E{z 

K^V ^. . 1 + z dE{z)^ 



(54) 



where E{z) = H{z)/Ho. Note that we have changed the time derivative to the derivative with respect to the redshift 
z, by using the relation dt = —dz/[H{l + z)]. Integrating Eq. ([Ml) , the field (j) is known as a function of z. Inverting 
<p{z) to z((j)) and plugging it into Eq. (|55p. one can reconstruct the potential V with respect to (f> by using the 

information of the observationally known values of H{z) and H'{z) as well as VL^m ■ We caution, however, that the 
actual observational data (such as the luminosity distance) are obtained at discrete values of redshifts. Hence we need 
some smoothing process for reconstructing the potential V{(j)) and the field equation of state W(f,{z). This smoothing 
was already discussed in Sec. Ill Al 

From the viewpoint of particle physics, the quintessence energy density can be comparable to the background fluid 
density in the early Universe. It is possible to construct quintessence models in which the field energy density is 
proportional to the fluid density during radiation and matter eras. For the exponential potential (|48l) there is a fixed 
point giving a constant field density parameter il^ — 3(1 + u'm)/A^ with = wm, where wm is the fluid equation of 
state ^35j . This is called a scaling solution, which is stable for > 3(1 +wm)- During radiation and matter eras one 
has = 4/A^ and fl^ = 3/A^, respectively. The Big BangNucleosynthesis (BBN) places the bound fl^ < 0.045 (95 
% confidence level) around the temperature T = IMeV j21J|, which gives the bound A > 9.4. In this case, however, 
the scalar field does not exit to the accelerated fixed point given by Eq. 

There are several ways to allow a transition from the scaling regime to the epoch of cosmic acceleration. One of them 
is to introd uce a field pote ntial that beco mes shallow at late times, V{4>) = cie~'^^'^ + C2e~'^'^'^ with A^ > 3(1 + wm) 
and p? <2 [215j (see Refs. |2l4 1216| - |221 [ for related works). For this double exponential potential the field equation 
of state of the final attractor is given by w^j, — —\ + p^ 13. Another way is to consider multiple scalar fields with 
exponent ial p otentials, y(0i,^2) = [222l . l223j | . In this case the phenomenon called assisted 

infiation j224| occurs for the multi-field exponential potential, even if the individual field has too steep a potential to 
lead to cosmic acceleration. The scalar field equation of state finally approaches the value = — 1 + A^g/3, where 
Aeff = (1/A^ + 1/A|)~^/^ is smaller than each A^ [i = 1, 2). I n the presence of three assisting scalar fields it is possible 
to realize the observational bound 100(2: = 0) < —0.8 today l225ll. 

There is a class of models dubbed quintessential inflation [2261] in which a single scalar fleld is responsible for 
both inflation and dark energy. One example of quintessential inflation is given by V{(j)) = \{4>'^ + M^) for (f> < Q 
and V{(t>) = AM*/ [1 + (0/M)"] [n > 0) for > 0. In the regime (j> < with |^| > M the potential behaves as 
V{(f)) ~ X4>'^, which leads to inflation. In the regime (j> > with (j) ^ M one has V{(j)) ~ AM*"*""!^"", which leads the 
late-time cosmic acceleration. Since the potential does not have a minimum, the reheating after inflation proceeds 
with a gravitational particle production. Although this process is not efficient in gen eral, it may be possible to make 
the reheating more efficient under the instant preheating mechanism proposed in Ref. 227 1. See Refs. (228.-.239j about 
related works on quintessential inflation. 



B. k-essence 



Scalar fields with non-canonical kinetic terms often appear in particle physics. In general the action for such theories 
can be expressed as 



d^x. 
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2^ 



R + P((h,X] 



S 



M 



(56) 
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where P{(j>,X) is a function in terms of a scalar field (j) and its kinetic energy X = —{l/2)g^'^dfj_(j)di^<t>, and Sm is a 
matter action. Even in th e absence of the field potential V{(f>) it is possible to realize the cosmic acceleration due to 
the kinetic energy X |240l | . The application of these theories to dark energy was first carried out by Chiba et al. [s^ . 
In Ref. [ioj this was extended to more general cases and the models based on the action (1551) were named "k-essence" . 
The action ([56| includes a wide variety of theories listed below. 



(A) Low energy effective string theory 

Th e act i on o f low energy effective string theory in the presence of a higher-order derivative term (V(/))* is given 



by [mill 



I d^iV^ F{^)R + oomVc^f + a'S^Vc^f + 0{ol^) , (57) 



which is derived by the expansion in terms of the Regge slope parameter a' (this is related to the string mass 
scale Ms via the relation Ms = y^2/a'). The scalar field (f), dubbed a dilaton field, is coupled to the Ricci scalar 
R with the strength F{(f>). This frame is called the Jordan frame, in which the tilde is used in the action (|57p. 
Under a conformal transformation, y^j, — F ((f))g ^^i, , we obtain the action in the Einstein frame [240j | 



1 

2k 



where K{<t)) = i{F^^lFf - 2uj/F and £(0) = 2a' B{4>) / . 



^R + K{(l))X + L{(t))X' 



(58) 



(B) Ghost condensate 

The theories with a negative kinetic energy —X generally suffers from the vacuum instability j243l |244| . but 
the presence of the quadratic term X^ can evade th is pr oblem. The model constructed in this vein is the ghost 
condensate model characterized by the Lagrangian '24511 

P = -X + XVM^ , (59) 

where M is a constant. A more general version of this model, called the dilatonic ghost condensate |246l |. is 

P = -X + e^'^'^'X^/M* , (60) 

which is motivated by a dilatonic higher-order correction to the tree-level action [as we have discussed in the 
case (A)]. 

(C) Tachyon 

A tachyon field appears as an unstable mode of D-branes [non-Bog omornyi-Prasad-Sommerfield (non-BPS) 
branes]. The effective 4-dimensional Lagrangian is given by [24714249 1 



P= -Vi(j))\/l-2X , (61) 

where V{(j)) is a potential of the tachyon field (f>. While the tachyon model is difficu lt to be compatible with 
infiation in the early Universe because of the problem for ending i nflation j250l - [25^ . one can use it for dark 
energy provided that the potential is shallower than V{(j)) = Vo(j)~'^ |253l - l258| . 

(D) Dirac-Born-Infcld (DBI) theory 

In the so-called "D-cceleration" mechanism in which a scalar field <j> parametrizes a direction on the approximate 
Coulomb branch of the system in A/" = 4 supersymmetric Yang Mills theory, the field dynamics can be described 
by the D BI action for a probe D3-brane moving in a radial direction of the Anti de Sitter (AdS) space-time 
[259l l26d |. The Lagrangian density with the field potential V{(j)) is given by 



P = Vl - 2/(0)A + /(0)-i - U(^) , (62) 

where /((/>) is a warped factor of the AdS throat. In this theory one can realize the acceleration of the Universe 
even in the r egime wh ere 2f{(j))X is close to 1. The application of this theory to dark energy has been carried 
out in Refs. f26ll - l263j . 
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For the theories with the action ([56)) . the pressure and the energy density of the field are = P and 
Pcf, = 2XP X ~ P, respectively. The equation of state of k-essence is given by 

As long as the condition |2XPjs:| ^ |P| is satisfied, can be close to —1. In the ghost condensate model ((59)) we 
have 

l-X/M* 
= 1 - ' 

which gives —1 < < —1/3 for 1/2 < X/M^ < 2/3. In particular the de Sitter solution {wtf, = —1) is realized at 
X/M'^ = 1/2. Since the field energy density is — M^/A at the de Sitter point, it is possible to explain the cosmic 
acceleration today for M ~ 10~^ eV. 

In order to discuss stability conditions of k-essence in the Ultra- Violet (UV) regime, we decompose the field into the 
homogenous and perturbed parts as (j){t, x) — (j)o{t) + 5(j){t, x) in the Minkowski backgro und a nd derive the Lagrangian 
and the Hamiltonian for perturbations. The resulting second-order Hamiltonian reads [246j | 

on = (P^x + 2XP^xx) — hPx — ^ P00 ^ ■ (65) 

The term P^^^ is related with the effective mass of the field, which is unimportant in the UV regime as long as the 
field is responsible for dark energy. The positivity of the first two terms in Eq. (|65p leads to the following stability 
conditions 

Px + 2XPxx > , Px >0. (66) 

The phantom model with a negative kinetic energy —X with a potential V{(j)), i.e. P = —X — V{(l)), do not satisfy 
the above conditions. Although the phantom model with P = —X — V{(j)) can lead to the background cosmological 
dynamics allowed by SN la observations {w^ < — 1) (243l 1264 1265| . it suffe rs fro m a catastr ophic particle production 
of ghosts and normal fields because of the instability of the vacuum i243l [24^ (see Refs. [iM |266- 280] for related 
works). This problem is overcome in the ghost condensate model ([SQ)) in which the conditions (|66|) are satisfied for 
X/M^ > 1/2 Thus the successful k-essence models need to be constructed to be consistent with the conditions 
while the field is responsible for dark energy under the condition |2XPx| ^ |P|. 

The propagation speed Cs of the field is given by |282j | 

2 P-P.X P,X 

which is positive under the conditions (1661) . The speed Cs remains sub- luminal provided that 

Pxx>0. (68) 



This condition is ensured for the models ([59| . (|60)) . (|6T|) . and (|62|) . 

There are some k-essence models proposed to solve the coincidence problem of dark energy. One example is [10, |4l| 

P = ^ (-2.01 + 2^/1 + X + 3 • 10""X^ - IQ-'^'^X^^ . (69) 

In these models the solutions can finally approach the accelerating phase even if they start from relatively large values 
of the k-essence energy density 51^ in the radiation era. In such cases, howe ver, there is a period in which the sound 
speed becomes superluminal before reaching the accelerated attractor j283l |. Moreover it was shown that the basins 
of attraction of a radiation scaling solution in such models are restricted to be very small j284| . We stress that these 
problems arise only for the k-essence models constructed to solve the coincidence problem. 



^ It is possible that the dilatonic ghost condensate model crosses the cosmological constant boundary = —1, but the quantum instability 
probl em is present in the region < — 1. This crossing does not occur for the single-field k-essence Lagrangian with a linear function 
oiX [28lll . 
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C. Coupled dark energy 

Since the energy density of dark energy is the same order as that of dark matter in the present Universe, this imphes 
that dark energy may have some relation with dark matter. In this section we discuss the cosmological viabihty of 
coupled dark energy models and related topics such as scaling solutions, the chameleon mechanism, and varying a. 

1. The coupling between dark energy and dark matter 

In the flat FLRW cosmological background, a general coupling between dark energy (with energy density pde and 
equation of state wde) and dark matter (with energy density pm) may be described by the following equations 

PDE + 3i?(l + u;DE)pDE - -P, (70) 

Pm + SHpm = +/3, (71) 

where /3 is the rate of the energy exchange in the dark sector. 

There are several forms of couplings proposed so far. Two simple examples are given by 

(A) I3 = kQp^^, (72) 

(B) p = aHp^, (73) 

where Q and a are dimensionless constants . The coupling (A) arises in scalar-tensor theori es after th e conformal trans- 
formation to the Einstein frame [285l - l288i . In general the coupling Q is fi eld-dependent |289l l29Clj | , but Brans-Dicke 
theory [29 ij (including f{R) gravity) gives rise to a constant coupling [295]. The coupling (B) is more phenomenolog- 
ical, but this form is useful to place observational bou nds from the cosmic exp ansion hi story. Several a uthors st udied 
other couplings of the forms (3 = {a„^pm + aBEPDE)H (293l - [295t . /3 = a^DE-ff (296l429i |. and /3 = Tpm p99l - l30Ti . See 
also Refs. I302l-l318ll for related works. 



(A) The coupling (A) 

Let us consider the coupling (A) in the presence of a coupled quintessence field with the exponential potential (|48|) . 
We assume that the coupling Q is constant. Taking into account radiation uncoupled to dark energy {pr oc a^*^) the 
Friedmann equation is given by 3H^ — k^{pbe + Pm + Pr): where pde = 0^/2 + 1^(0)- Introducing the dimensionless 
variables xi = k4'/{\/6H), X2 — kVV /{VSH), and X3 = Ky/p^/{\/3H) as in Sec. IIV Al we obtain 

x[ = -3x1 + ^^^l + ^^1 {i^'ixl~3xl+xl)-^Q{l-xl~xl-xl), (74) 

and the same differential equations for X2 and 2:3 as given in Eqs. (1461) and (|47p . For this dynamical system there 
is a scalar-field dominated fixed point given in Eq. (j49p as well as the radiation point (xi,a;2,a;3) = (0,0, 1). I n the 
presence of the coupling Q, the standard matter era is replaced by a "(/)-matter-dominated epoch ((/)MDE)" j287l | 
characterized by 

(xi,a;2,X3) = (-y6Q/3,0,0), fi^ = 2QV3 , = \ . (75) 

Defining the effective equation of state 

u-eff = -1 - 2H/{3H^) , (76) 

one has WcS = '2,Q^ /3 for the 0MDE, which is different from in the uncoupled case. Provided that 2(5^/3 < 1, the 
(/)MDE is a saddle followed by the accelerated point (|49|) j287| . 

The evolution of the scale factor during the (/iMDE is given by a oc t^/i^+'^Q ) ^ which is different from that in 
the uncoupled quintessence. This leads to a change to the CMB shift parameter define d in Eq. ([2^ . From the 
CMB likelihood analysis, the strength of the coupling is constrained to be IQj < 0.1 [287| . The evolution of matter 
density perturbations is also subject to change by the effect of the coupling. Under a quasi-static approximation on 
sub-horizon scales the matter perturbation 5m obeys the following equation [319, 320] 

5 m + (2i? + Q(i>)5m — 47rG' 

off Pm^m 

~0, (77) 
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where the effective gravitational coupling is given by GcS = (1 + '2Q^)G. During the (/)MDE one can obtain the 
analytic solution to Eq. ([77]) . as Sm oc _ Hence the presence of the coupling Q leads to a larger g rowt h rate 



relative to the uncoupled quintessence. We can parameterize the growth rate of matter perturbations, as |321l 



H6„ 



(78) 



where ilm = n^Pm/i'^H^) is the density para meter of n on-relativistic matter. In the ACDM model the growth index 
7 can be approximately given by 7 ~ 0.55 [322l l323j |. In the coupled quintessence the growth rate can be fitted 
to the numerical solution by the formula / = (0^)^(1 + cQ^), where c = 2.1 and 7 = 0.56 are the best-fit values 
^24]. Using the galaxy and Lyman-a power spectra, the growth index 7 and the coupling Q are constrained to be 
7 = 0.6 and \Q\ < 0.52 (95 % confidence level), respectively. This is weaker than the bound coming from the CMB 
constraint |324| . We also note that the equatio n for matter perturbations has been derived for the coupled k-essence 
scenario with a field-dependent coupling Q{(f>) [325l l326l |. In principle it i s pos sible to reconstruct the coupling from 
observations if the evolution of matter perturbations is known accurately [3271 . 

There is another interesting coupled dark energy scenario called mass varying neutrino [3281433 j | in which the mass 
nil, of the neutrino depends on the quintessence field (j). The energy density pi, and the pressure of neutrinos can 
be determined by assuming a Fermi-Dira c dist ribution with the neglect of the chemical potential. It then follows that 
the field cj) obeys the equation of motion [33 1| 



4> + + ^0 = -kQ((/))(p^ - 3P^) , 



(79) 



where Q{(j)) — fJilL^iW) ^ jj^ i^j^g relativistic regime in which the neutrino mass rrii, is much smaller than the neutrino 
temperature T^, the r.h.s. of Eq. (|79p is suppressed because of the relation p^, ~ 3Py. In the non-relativistic regime 
with rrii, 3> the pressure P^, is much smaller than the energy density p^, in which case the field equation (j79p mimics 
Eq. (I7n|) with the coupling (A). In the mass varying neutrino scenario the field-dependent mass of neutrinos determines 
the strength of the coupling Q (</»). In the non-relativistic regime the neutrino energy density is approximately given 
by pi, ~ ni,m^{(l)), where rij, is the number density of neutrinos. Then the effective potential of the field is given 
by Vcs{(f>) — V{(j)) -\- n^m,^{(j))^ which gives rise to a minimum for a runaway quintessence potential V{(j)). Since the 
field equation of state in this regime is ~ —V{4))/[V{(/)) + ni^rrii, {(/))], one has ~ —1 for nj^rrii, <C V^(0). See 
Refs. [333l - l338t for a number of cosmological consequences of the mass varying neutrino scenario. 



(B) The coupling (B) 



Let us consider the coupling (B) given in Eq. (|73l) . For constant 5 one can integrate Eq. (|7ip as pm = p^m i^ + z)^ 
where z is a redshift. If the equation of state wde is constant, one obtains the following integrated solution for Eq. (j70p : 



3it;DE 



(1 + ^) 



3(1 + u)de) 



(1 + ^) 



3-a 



(80) 



The Friedmann equation, 37?^ 
E{z) = H{z)/Ho, as 



k^(pde 



Pm), gives the parametrization for the normalized Hubble parameter 



E^iz)^n'Slii + z) 



3(1+iode) 



1 - n 



(0) 
DE 



a + SwDE 



a(l-Kz)3(^+"°^)+3wDE(l + 2)' 



(81) 



This parametrization can be used to place observational constraints on the cou plin g a. The combined data analysis 
using the observational data of the 5-year Supernova Legacy Survey (SNLS) [IJl, the CMB shift parameter from 
the 3-year WMAP [20|, and the BAO [2^ shows that a and wde are constrained to be —0.08 < a < 0.03 and 
— 1.16 < wuE < —0.91 (95 % confidence level) |298j |. In Ref. ^300j it was shown that, for constant wbe, cosmological 
perturbations are subject to non-adiabatic instabilities in the early radiation era. This pr oble m can be alleviated by 
considering the time-dependent wbe satisfying the condition wde > —4/5 at early times j317j |. 

It is also possible to extend the analysis to the case in which the coupling a varies in time. Dalai et al. |296l | 



assumed the scaling relation, p-DE/Pm — [p 
expressed in the form 



de/a 



a^, where ^ is a constant. For constant weie the coupling a is 



a{z) 



^^^l + {i-n^^i){i + z) 



(0). 



(82) 
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where ao = —{C + 3wDE)i^DE- If ^ > 0, a{z) decreases for higher z. In this case the Hubble parameter can be 
parametrized as 

E'iz) = (1 + z f [l ~ d^^^ + + z)-^\ . (83) 

The combined data analysis using the SNLS, WMAP 3-year, and the BAO gives the bounds —0.4 < ao < 0.1 and 
-1.18 < tWDE < -0.91 (95 % confidence level) j298j |. The ACDM model (a = 0, wde = -1) remains a good fit to the 
data. 

2. Coupled dark energy and coincidence problem 

In the coupled quintessence with an exponential potential {V{4>) = Vqc"'''^'^ with A > 0) the </)MDE scaling solution 
wit h ^<j > = 2 Q^/3 = constant replaces the standard matter era. In this model there is another scaling solution given 
by [287[l306i 



f ^ ( ^ 2Q(Q + A) + 3 \ Q{Q + X) + 3 

^"^'"^'"^^^1,^(QTa)'\/ 2iQ + xr ''J ' ""'^ {Q + Xr ' 

Q{Q + X) Q 
"^ = "Q(Q + A) + 3' "'=« = -oTa- 

In the presence of the coupling Q the condition for the late-time cosmic acceleration, < —1/3, is satisfied for 
Q > A/2 > or Q < —A < 0. Then the scaling solution ([M|) can give rise to the global attractor with 51^ ~ 0.7. 
The (j)MDE solution (1751) followed by the accelerated scaling solution ([M)) may be used for alleviating the coincidence 
prob lem, because dark energy and dark matter follow the same scaling relation from the end of the radiation era. In 
Ref. 287] it was shown, however, that the coupled quintessence with an exponential potential does not allow for such 
cosmological evolution for the viable parameter space in the {Q, A) plane consistent with observational constraints. 

It is possible to extend the analysis to coupled k-essence models described by the action (f56| in the presence of the 
coupling ()72p . From the requirement that the density parameter (7^ 0) and the equation of state w^, are constants 
to realize scaling solutions, one can show that the Lagrangian density takes the following form [246i . i339i] 

P^XgiXe''^'''), (85) 

where A is a constant, and g is an arbitrary function in terms oiY = Xe'^^'^. The result is valid for constant Q, 
but it can be generalized to a field-dependent coupling Q{(j)) t34(j]. The quintessence with an exponential potential 
{P = X — ce""^*^) corresponds to g{Y) — 1 — c/Y. Since the dilatonic ghost condensation model corresponds to 
g{Y) = — 1 + F/M^, this model has a scaling solution. We can also show that the tach yon field with the Lagrangian 
density P = —V{(l))y/1 — 2X also has a scaling solution for the potential V{(j)) oc (/)~^ |257| . Even when the Hubble 
parameter squared is proportional to the energy density p", it is possible to obtain the Lagrangian density having 
scaling solutions in the form P = X-^/"-n(Xe ^'^^'^) [339]. The cosmological dyna niics of sc aling solutions in such cases 
(including the high - energy r egime (34ll l342j in the Randall-Sundrum scenario |343l . 134^ ) have been discussed by a 
number of authors p5l - [347j . 

For the Lagrangian density (|55|) there are two fixed points relevant to dark energy. Defining the dimensionless 
variables x = (j)/{\/QH) and y = e"^'^^^ / {^/SH) (in the unit of = 1), they are given by j348j | 

• (A) Scalar-field dominated solution 

A A^ 
^A = —7^-—, fl0 = l, wcs = w^, = -1 + —— . (86) 



(B) Scaling solution 



VQ „ QiQ + X) + Wx 



2 



2(Q + A)' " (Q + A) 

Q Q{Q + X) 

"'=^ = -QTa' "^ = " Q(Q + A) + 3Px - ^''^ 
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The points (A) and (B) are responsible for the cosmic acceleration for (A) \^ jP.x < 2 and (B) Q > A/2 > or 
Q < —A < 0, respectively. From the stability an a,lysis about the fixed points it follows that, when the point (A) is 
stable, the point (B) is not stable, and vice versa |348j |. 

The (pMDE corresponds to a fixed point at which the kinetic energy of the field dominates over the potential energy, 
i.e. X ^ and y — 0- Since the quantity Y = Xe^^'^ can be expressed as y = x^ /y^, the function g{Y) cannot be 
singular at y = for the existence of the (/)MDE. Then the function g{Y) should be expanded in negative powers of 
y, i.e. 

g{Y) = CO + ^ c„r-" - CO + ^ c^ivV^^T , (88) 

n>0 n>0 

which includes the quintessence with an exponential potential. For this form of g(Y), there is the following 0MDE 
point (C): 

[xc,yc)=\ — 5 — ,0 , Hp^Wcft^- — , = (89) 
\^ 3co J 3co 

together with the purely kinetic point {x, y) = (±l/.yco, 0) and 17^ = 1 for co > 0. 

An ideal cosmological trajectory that alleviates the coinci denc e problem of dark energy should be the 0MDE (C) 
followed by the point (B). However, it was shown in Ref. |340l | that such a trajectory is not allowed because the 
solutions cannot cross the singularity at a; = as well as another singularity associated with the sound speed. For 
example, when cq > 0, we find that > and < for Q > A/2 > 0, whereas < and > for Q < —A < 0. 
These points are separated between the line a; = at which the function ([55)1 diverges. The 0MDE solution chooses 
the accelerated point (A) as a final attractor. The above discussion shows that the coincidence problem is difficult to 
be solved even for the general Lagrangian density ()85|) that has scaling solutions. This problem mainly comes from 
the fact that a large coupling Q required for the existence of a viable scaling solution (B) is not compatible with 
a small co uplin g Q required for the existence of the 0MDE. We need a rapidly growing coupling to realize such a 
transition 134911. 



3. Chameleon mechanism 



If a scalar field (j) is coupled to baryons as well as dark matter, this gives rise to a fifth force interaction that can 
be constrained experimentally. A large coupling of the order of unity arises in modified gravity theories as well as 
superstring theories. In such cases we need to suppress such a strong interaction with baryons for the compatibility 
with local gravity experiments. There is a way to screen the fifth force under the so-called chameleon mechanism 
[350ll35i| in which the field mass is different depending on the matter density in the surrounding environment. If the 
field is sufficiently heavy in the regions of high density, a spherically symmetric body can have a "thin-shell" around 
its surface such that the effective coupling between the field and matter is suppressed outside the body. 

The action of a chameleon scalar field with a potential V{(t>) is given by 

+ j d\CM{g\:lM[l), (90) 

where g is the determinant of the metric g^u (in the Einstein frame) and Cm is a matter Lagrangian with matter 

fields ^^'j coupled to a metric g'lil. The metric is related to the Einstein frame metric g^^, via g'lil — e^'^'^^'^g^^, 
where Qi are the strengths of the couplings for each matter component with the field (p. The typical field potential 
is chosen to be of the runaway type (such as ^(0) = Af''+"(/)""). We also restrict the form of the potential such that 
-> oo as (/i -> 0. 

Varying the action ([90)) with respect to (f), we obtain the field equation 

- = - E '^Q.e^'^Q'^sf.^rW , (91) 

i 

where T^^} = — (2/ v^— g^'^ )5Cm /Sg^^ is the stress-energy tensor for the i-th form of matter. For non-relativistic matter 

we have g'^^^T^^} — —pi, where pi is an energy density. It is convenient to introduce the energy density pi = piC^"'^''^, 
which is conserved in the Einstein frame. In the following, let us consider the case in which the couplings Qi are the 



d^x 



J-^R^^gt^-d^<l>d,<l>-V{<l>) 
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same for all species, i.e. Qi = Q. In a spherically symmetric space-time under the weak gravitational background 
(i.e. neglecting the backreaction of gravitational potentials), Eq. ([?!]) reads 

d20^2d0^dK£^ 

dr'^ r dr d(f) ' 

where r is a distance from the center of symmetry, and Vcn is the effective potential given by 

Kff(0) = l^(0)+e''«*p, (93) 

and p = Pi. For the runaway potential with < the positive coupling Q leads to a minimum of the effective 
potential. In f[R) gravity the negative coupling [Q — — 1/-\/6) gives rise to a minimum for the potential with V^^ > 
(as we wiU see in Sec. IV A3|) . 

We assume that a spherically symmetric body has a constant density p — pA inside the body {r < Tc) and that the 
energy density outside the body (r > Tc) \s p = pb- The mass of the body and the gravitational potential $c at 
the radius Tc are given by Mc = (47r/3)rj^pA and $c = GMc/vc, respectively. The effective potential Vcs{(t)) has two 
minima at the field values (f)A and satisfying V^g{4'A) = and V^'fj((^_B) = 0, respectively. The former corresponds 
to the region with a high density that gives rise to a heavy mass squared m\ = V^g{(j)A), whereas the latter to the 
lower density region with a lighter mass squared = ^/^((/ib). When we consider the "dynamics" of the field cj) 
according to Eq. (I92p we need to consider the inverted effective potential {—Vcs) having two maxima at (f) = (pA and 

The boundary conditions for the field are given by ^(r = 0) = and (j){r oo) = (pB- The field is at rest at 
r — and begins to roll down the potential when the matter-coupling term kQpas'^'^^ becomes important at a radius 
ri in Eq. (|92|). As long as ri is close to Vc so that Arc = ~ ri ^ rc, the body has a thin-shell inside the body. 
Since the field acquires a sufficient kinetic energy in the thin-shell regime (ri < r < rc), it climbs up the potential 
hill outside the body (r > r^). The field profile can be obtained by matching the solutions of Eq. (1^^ at the rad ius 
r — ri and r — rc- Neglecting the mass term tub, we obtain the thin-shell field profile outside the body [350l4352{ 

(/)(r) ~ 0s , (94) 

K r 

where 

n - k(0b - Pa) , . 

Qcs = Sgeth , Eth = ■ (95) 

Here eth is called the thin-shell parameter. Under the conditions lS.rc/rc ^ 1 and l/(m^rc) ^ 1, the thin-shell 
parameter is approximately given by eth — At^/ + 1/ {mArc) 352] . As long as eth ^ 1 the amplitude of the effective 
coupling Qcff can be much smaller than 1. Hence it is possible for the large coupling models {\Q\ = 0(1)) to be 
consistent with local gravity experiments if the body has a thin-shell. 

Let us study the constraint on the thin-shell parameter from the possible violation of the equivalence principle (EP). 
The tightest bound com es fro m the solar system tests of weak EP using the free-fall acceleration of Moon (aMoon) and 
Earth (a®) toward Sun [35 1| . The experimental bound on the difference of two accelerations is given by [353i . i354i] 

IflMoon^Q-©! ,-,„_i3 



(flMoon + ae)/2 



< 10-'^ (96) 



If Earth, Sun, and Moon have thin-shells, the field profiles outside the bodies are given by Eq. ([94]) with the replacement 
of corresponding quantities. The acceleration induced by a fifth force with the field profile (/)(r) and the effective 
coupling Qefi is a*^^"^ = |(5effV^(r)|. Using the thin-shell parameter eth,® for Earth, the accelerations a® and flMoon 
toward Sun (mass Mq) are [35 1| 



GMo 



© 

2,2 



l + 18Q2e?h,eg 



1 + 18Q^et\,e 



©y^Moon 



(97) 
(98) 



where ~ 2.1 x 10 ^, $0 ~ 7.0 x 10 and $Moon — 3.1 x 10 are the gravitational potentials of Sun, Earth 
and Moon, respectively. Hence the condition (IMl) translates into 

eth,® < 8.8 X 10-VlQI . (99) 
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Since the condition \(j)B \ ^ I'f'Al is satisfied for the field potentials under consideration, one has eth,( 
from Eq. (|95p . Then the condition ([M)) translates into 



|K(/.B,e| < 3.7 X 10- 



-15 



k0b/(6Q$®) 



(100) 



For example, let us consider the inverse power-law potential V{(p) — M"'+"0 In this case we have (/is,® — 
[(n/(3)(Mj^j/pB)(M/Mpi)"+'*]i/("+i)Mpi, where we recovered the reduced Planck mass Afpi = 1/k. For n and Q 

of the order of unity, the constraint (jlOOl) gives M < lO^(^^"+^-^'')/("+'')Mpi. When n — 1, for example, one has 
M < 10~^ eV. If the same potential is responsible for dark energy, the mass M is constrained to be larger than this 
value [Hi]. For the potential V{(j)) = M'' exp(M"/0"), however, we have that ^ + M'*+"0-" for (/) > M, 

which is respo nsible for d ark energy fo r M » 10~'^ eV . This can be compatible with the mass scale M constrained 



by (Unni) [3555. See Refs. 
chameleon field. 



359l l359l - l363l | for a number of cosmological and experimental aspects of the 



4- Varying a 

So far we have discussed the coupling between dark energy and non-relativistic matter. In this section we discuss 
the case in which dark energy is coupled to an electromagnetic field. In fact, a temporal variation of the effective 
fine structure "constant" a has been reported by a number of authors. The variation of a constrained by the Oklo 
natural fission reactor is given b y —0 .9 x 10^^ < Aa/a = {a — ao)/ao < 1.2 x 10~^ at the redshift z sa 0.16, 
where ao is the value of a today |364| . From the absorption line spect ra o f dist ant quasars we have the constraints 
Aa/a = (-0.574 ±0.102) x 10"^ over the redshift range 0.2 < z < 3.7 [sel 13665 and Aa/a = (-0.06 ±0.06) x 10"^ 
for 0.4 < z < 2.3 [367.] . Although the possibility of systematic errors still remains [368], this may provide important 
implications for the existence of a light scalar field related with dark energy. 

The Lagrangian density describing such a coupling between the field and the electromagnetic field Ffj^i, is given 

by 

£f{^)^-^Bf{cI>)F^.F^-'. (101) 

The coupling of the form Bp{(j)) — e~'»''("^~"^") was originally introduced by Bekenstein 1369*], wh ere C is a coupling 
constant and (j)a is the field value today. There are also other choices of the coupling, see e.g., Refs. (STOmTTt . The fine 



structure "constant" a is inversely proportional to Bp{(l)), so that this can be expressed as a = ao/ Bp{4>), where ao 
is the present value. The exponential coupling Bf{4>) = e~'''^^^~^°'^ has a linear dependence Bf{4>) ~ 1 — CK{(f> — <po) 
in the regime \Ck,{4> — '/'o)| ^ 1, so that the variation of a is given by 

^ = ^-C-(0-0o). (102) 
a ao 



Using the constraint Aa/a ~ — 10 ^ around 2 = 3 j366l] obtained from quasar absorption lines, the coupling ^ can be 
expressed as 

10"^ 

C^--7 TT. -i ^- (103) 

K0(Z = 3) — K(p{Z = 0) 

Let us consider the case in which the scalar field has a power-law dependence in terms of the scale factor a, i.e. 

<i>^<PoaP ^U^ + z)-P . (104) 

In fact, in the so-called tracking regime of the matter-dominated era [11], the inverse power-law potential V{4>) = 
j^i+n^-n gjygg j.jgg ^ constant field equation of state: = —2/{n + 2), which corresponds to the field evolution 
([TW)) with p = 3/(n + 2). Using Eq. pTO)) . the coupling C in Eq. (IT^ reads 

1 0^^ 

(:^——{nc^,y\ (105) 

Since k0o is of the order of unity in order to realize the present cosmic acceleration (ssj , the coupling C, is constrained 
to be C ~ 10~^ for p of the order of 1. 

The above discussion is valid for the potentials having t he_s olution ()104p in the tracking regime. The variation of 
a for other quintessence potentials was discussed in Ref. |3 76j . There is also a k-essence model in which a tachyon 
field is coupled to electromagnetic fields j378|. 
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D. Unified models of darlc energy and darlc matter 



There are a number of works to explain the origin of dark energy and dark matter using a single fluid or a single 
scalar field. Let us first discuss the generalized Chaplygin gas (GCG) model as an example of a single fluid model 
[3 ■ this model the pressure P of the perfect fluid is related to its energy density p via 



P 



-Ap- 



(106) 



where A is a positive constant. The original Chaplygin gas model corresponds to a = 1 [44|. 

Plugging the relation (jl06|) into the continuity equation p + 3H{p + F) = 0, we obtain the integrated solution 



A 



B 



a3(l+a) 



(107) 



where B is an integration constant. In the early epoch (a ^ 1) the energy density evolves as p cx a~^, which means 
that the fluid behaves as dark matter. In the late epoch (a ^ 1) the energy density approaches a constant value 
^i/(a+a) hence the fluid behaves as dark energy. A fluid with the generalized Chaplygin gas therefore interpolates 
between dark matter and dark energy. 

Although this model is attractive to provide unified description of two dark components, it is severely constrained 
by the matter power spectrum in large-scale structure. The gauge-invariant matter perturbation 5m with a comoving 
wavenumber k obeys the following equation of motion [46i] 



(2 + 3c2 - 6w) H5m - 



Ih\1 - 6c2 - 2,w^ + 8w)-('^ 
2 \ a 



5„,, = . 



where w = P/ p is the fluid equation of state, and Cg is the sound speed given by 



2 dP 
c„ = — — = ~aw . 
dp 



(108) 



(109) 



Since w — ^ and — >■ in the limit z ^ 1, the sound speed is much smaller than unity in the deep matter era and 
starts to grow around the end of it. Since w is negative, cj. is positive for a > and negative for a < 0. 
From Eq. (jlOSp the perturbations satisfying the following condition grow via the gravitational instability 



3 faH 
<2 - 



(110) 



When |c^| > (3/2) [aH/k) , the perturbations exhibit either rapid growth or damped oscillations depending on the 
sign of c^. The violation of the condition (IllOp mainly occurs around the present epoch in which \'w\ is of the order 
of unity and hence \c1\ ~ \a\. The smallest scale relevant to the galaxy matter power spectrum in the linear regime 
corresponds to the wavenumber around k = 0.1 /i Mpc^^. Then the constraint ()110|) gives the upper bound on the 
values of lal [4611: 



< 10 



-5 



(111) 



Hence the generalized Chaplygin gas model is hardly distinguishable from the ACDM model. In particular the original 
Chaplygin gas model (a = 1) is excluded from the observa tion s of la rge-scale structure. Although nonlinear clustering 
may change the evolution of perturbations in this model j379L l380l | , it is unlikely that the constraint (|111|) is relaxed 
signiflcantly. 

The above conclusion comes from the fact that in the Chaplygin gas model the sound speed is too large to match 
with observations. There is a way to avoid this problem by adding a non-adiabatic contribution to Eq. (|108|) to make 
Cs vanish [381J. It is also possible to construct unified models of dark energy and dark matter using a purely kinetic 
scalar field |382l |. L et us consider k-essence models in which the Lagrangian density P{X) has an extremum at some 
value X = Xn, e 



382 [ 



P = P^ + P2{X - Xaf 



(112) 



The pressure P4, — P and the energy density = 2XP^x — P satisfy the continuity equation p0 -t- 3H{p^ + P^) = 0, 
i.e. 



(Fx + 2XPxx) X + GHP^xX = 



(113) 
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The solution around X = Xq can be derived by introducing a small parameter e = {X — Xq)/Xq. Plugging 
Eq. (11121) into Eq. (|113l) . we find that e satisfies the equation e = ~3He at linear order. Hence we obtain the 
solution X ~ Xq [l + ei(a/ai)^'^] , where ei and ai are constants. The validity of the above approximation demands 
that ei{a/ai)~^ ^ 1. Since ~ Pq and ~ — Pq + 4P2-'fo^i(a/ai)^'^ in the regime where X is close to Xq, the 
field equation of state is given by 



4^2 9 



(114) 



Since — >■ — 1 at late times it is possible to give rise to the cosmic acceleration. One can also realize « during 
the matter era, provided that the condition 4P2-^o/l^ol ^ 1 is satisfied. The sound speed squared defined in Eq. ([57)) 
is approximately given by 



1 



(115) 



which is much smaller than unity. Hence the large sound speed problem can be evaded in the model (jll2p . In 
Rcf. 383] it was shown that the above purely k-essence model is equivalent to a fluid with a closed-form barotropic 
equation of state plus a constant term that works as a cosmological constant to all orders in structure formation. See 
Refs. [384 - .388,] for generalized versions of the above model. 



V. MODIFIED GRAVITY MODELS 



There is another class of dark energy models in which gravity is modifled from General Relativity (GR). We review 
a number of cosmological and gravitational aspects of f{R) gravity, Gauss-Bonnet gravity, scalar-tensor theories, and 
a braneworld model. We also discuss observational signatures of those models to distinguish them from other dark 
energy models. 



A. f{R) gravity 

The simplest modification to GR is f{R) gravity with the action 

^ " 2^/ ^^'•^^/(^) + / d^^CMi9^.,'^M) , (116) 

where / is a function of the Ricci scalar R and Cm is a matter Lagrangian for perfect fiuids. The Lagrangian Cm 
depends on the metric and the matter field s W e do not consider a direct coupling between the Ricci scalar 
and matter (such as fi{R)CM studied in Refs. |389l - l391] ). 



1. Viable f{R) dark energy models 

In the standard variational approach called the metric formalism, the affine connection s T^^ are related with the 
metric 5^ 

u [QTI ] ■ In this formalism, the field equation can be derived by varying the action (|116|) with respect to 

F{R)R^,{g) - ^f{R)g^, - V^V,P(P) + gf^^DFiR) = , (117) 

where F{R) = df /dR, and T^jy = — (2 / y/—g)SCM / Sg^'^ is the energy-momentum tensor of matter. Note that there 
is another way for the variation of the action called the Palatini formalism in which the metric and the connections 
are treated as independent variables. In Sec. IV A 41 we shall briefly mention the application of Palatini f{R) gravity 
to dark energy. The trace of Eq. (I117P is given by 

3DF{R) + F{R)R-2f{R)^ K^T, (118) 

where T = y^T^jjy = — pm + 3Pm- Here pM and Pm are the energy density and the pressure of matter, respectively. 
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The de Sitter point corresponds to a vacuum sohition at which the Ricci scalar is constant. Since DF{R) = at 
this point, we obtain 

- 2/(i?) = . (119) 

The model f{R) — aR^ satisfies this condition and hence it gives rise to an exact de Sitter solution. In fact the first 
model of infiation proposed by Starobinsky [113] corresponds to f{R) = R + aR^, in which the cosmic acceleration 
ends when the term aR^ becomes smaller than R. Dark energy models based on f{R) theories can be also constructed 
to realize the late-time de Sitter solution satisfying the condition (I119p . 

The possibility of the late-time cosmic acceleration in f{R) gravity was first suggested by Capozziello 14711 in 2002. 
An f(R) dark energy model of the form /(i?) = R — (n > 0) was proposed in Refs. [49l - l52l | (see also 

Refs . [392r,397j ). but it became clear that this model suffers from a number of problems such as the m atter instability 
[398| . absence of the matter era |399l |40C| . and inability to satisfy local gravity constraints [40114406 [. This problem 
arises from the fact that f^RR < in this model. 

In order to see why the models with negative values of f^RR, are excluded, let us consider local fiuctuations on a 
background characterized by a curvature Rq and a density po- We expand Eq. (|118p in powers of fluctuations under 
a weak field approximation. We decompose the quantities F(R), 17^1,, and T^^ into the background part and the 

perturbed part: R — i?'"' + SR, F = F^^\l + dp), g^,, — r/^i, -I- h^^, and T^,, — xjS) + 5T^^, where we have used 

the ap proximat ion that g'l^} corresponds the metric rj^i, in the Minkowski space-time. Then the trace equation (jllSp 
reads Eol, Eol 



Sf + A/ Sf = 



3F(") 



<5r, 



where ST = tj'^'ST^^, and 



i(i?(o)) 



(120) 



(121) 



Here the quantity m = Rf,RRjf,R characterizes the deviation from the ACDM model {f{R) = R — 2A). In the 
homogeneous and isotropic cosmological background (without a Hubble friction), (5^^^ is a function of the cosmic time 
t only and Eq. (|120l) reduces to 



5f 



5f = 



3F(o) 



(122) 



where p = —ST. For the models where the deviation from the ACDM model is small, we have m(i?*^°^) <C 1 
so that |M^| is much larger than R'-^K If A/^ < 0, the perturbation Sf exhibits a violent instability. Then the 
condition ~ /,fl(^^°^)/(3/,i?fl(i?^°^)) > is needed for the stability of cosmological perturbations. We also require 
that /,_r(-R^°-') > to avoid anti-gravity (i.e. to avoid that the graviton becomes a ghost). Hence the condition 
/,flfl(i?'"') > needs to hold for avoiding a tachyonic instability associated with the negative mass squared [4071 - 141 1| . 

For the consistency with local gravity constraints in solar system, the function f{R) needs to be close to that in the 
ACDM model in the region of high density (in the region where the Ricci scalar R is much larger than the cosmological 
Ricci scalar Rq today). We also require the existence of a stable late-time de Sitter point given in Eq. (|119p. From 
the stability analysis about the de Sitter point, one can show that it is stable for < m = RJ^rrJ/m < 1 412| - |414 |. 
Then we can summarize the conditions for the viability of f{R) dark energy models: 

• (i) /m > for i? > Ro. 

• (ii) f,RB. > for R>Rn. 

• (in) f{R) R-2A for R > Rq. 

• (iv) < R/mr/ f.R < 1 at the de Sitter point satisfying Rf_R — 2/. 

The examples of viable models satisfying all these requirements are 415l - l417 | 

{R/Rc?" 



(A) f{R) = R- /ii?. 



[R/Rcf'' + 1 

2 / d2' 



with n, p,, Rc > , 



(B) f{R) ^R- pR, 1 - (1 + RVR, 



with p, Rc > , 



(C) f{R) =R- pR^tanh {R/Rc) with p, Rc > , 



(123) 

(124) 
(125) 
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where ^, Rc, and n are constants. Models similar to (C) were proposed in Refs. |418l . l419l |. Note that Rc is roughly 
of the order of the present cosmological Ricci scalar Rq. If ^ i?c the models are close to the ACDM model 
(/(i?) ~ i? — fiRc), so that GR is recovered in the region of high density. The models (A) and (B) have the following 
asymptotic behavior 

/(i?)~i?-A*i?c[l-(i?/i?c)"'"] , (i?»i?c), (126) 

which rapidly approaches the ACDM model for n > 1. The model (C) shows an even faster decrease of m in the 
region R > R^. The model f{R) = R- fiRdR/RcY' (0 < n< 1) proposed in Refs. [iil [42^ is also viable, but it 
does not allow the rapid decrease of m in the region of high density required for the consistency with local gravity 
tests. 

For example, let us consider the model (B). The de Sitter point given by the condition (|119p satisfies 

^ 2[(l + a;2)»+i~l-(7i + l)a;2] ' ^ ' 

where xi = Ri/Rc and Ri is the Ricci scalar at the de Sitter point. The stability condition (0 < m < 1) at this point 
gives [416i] 

(1 + .t2)"+2 > 1 + (n + 2)xl + {n + l)(2n + 1)2;^ . (128) 

The condition (|128p gives the lower bound on the parameter fi. When n = 1 one has xi > and /i > Sa/S/S. Under 
Eq. (|128|) one can show that the conditions f^u > and /,_r_r > are also satisfied for R> Ri. 



2. Observational signatures of f{R) dark energy models 

In the fiat FLRW space-time we obtain the following equations of motion from Eqs. (|117p and (|118p : 

3Fij2 = K^pra + iFR- f)/2-SHF, (129) 
2FH = -K^prn - F + HF , (130) 

where, for the perfect fiuid, we have taken into account only the non-relativistic matter with energy density p,„. In 
order to confront f{R) dark energy models with SN la observations, we rewrite Eqs. (jl29p and (I130p as follows: 

3AH^ = (p™ + pue) , (131) 
-2AH - {p„, + pDE + Pde) , (132) 

where A is some constant and 

kVde = {1/2){FR- f)-iHF + ?,H''{A~F), (133) 
k^Pde = F + 2HF - {1/2){FR- f) - {3H^ + 2H){A~ F) . (134) 

By defining poE and Pde hi this way, one can easily show that the following continuity equation holds 

PDE + 3i?(pDE + PDE) =0. (135) 

We define the dark energy equation of state iude = Pd e/pd Ei which is directly related to the one used in SN la 



observations. From Eqs. (I13ip and (jl32l) it is given by [83|, 1421 



2AH + SAH^ w,s 
^°'=^ = -3Ai/2-«V™"l-(FM)17„' ^'''^ 

where = Pm / [iF H^) . The viable f{R) models approach the ACDM model in the past, i.e. — > 1 as i? — > oo. 
In order to reproduce the standard matter era in the high-redshift regime we can choose A = \ in Eqs. (I13ip and 
(|132l) . Another possible choice is A = Fq, where Fq is the present value of F. This cho ice is suitable if the deviation 
of Fq from 1 is small (as in scalar-tensor theory with a massless scalar field l422ll423|). I n both ca ses the equation 
of state wde can be smaller than —1 before reaching the de Sitter attractor [83l l415l |417| . l419l |424| . This originates 
from the fact that the presence of non-relativistic matter makes the denominator in Eq. (jl36p smaller than 1 (unlike 
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Refs. [425l l426l | in which the authors did not take into account the contribution of non-relativistic matter). Thus 
/(i?) dark energy models give rise to a phantom equation of state without violating stability conditions of the system. 
The mod els (A) and (B) are allowed from the SN la observations provided that n is larger than the order of unity 
[427l - l430t . 

The modification of gravity manifests itself in the effective gravitational coupling that appears i n the equ ation of 
cosmological perturbations. The full perturbation equations in f{R) gravity are presented in Refs. [43ll4433 [. When 
we confront /(i?) models with the observations of large-scale structure, the wavenumbers k of interest are sub-horizon 
modes with k/a 3> H. We can employ a so-called quasi-static approximation under which the dominant terms in 
perturbation equations correspond to those including fc^/a^, Sp m, and i9lU423l . l434| . Then the matter density 

perturbation 6m approximately satisfies the following equation [434l |435| | 



Sm + 2HSm - iTTGcSPmSm — , 

where pm is the energy density of non-relativistic matter, and 

G l + Amk^/{a^R) 



Gcs — 



f.Rl + 3mk^/{a^R) ' 



(137) 



(138) 



where m = Rf^RR/f.R- This approximation is accurate for viable f{R) dark energy models as long as an oscillating 
mode of the scalar-field degree of freedom is suppressed relative to the matter- induced mode |416l . l417i 1436 - 439]. 

In the regime where the deviation from the ACDM model is small such that mk^/{a^R) <^ 1, the effective gravita- 
tional coupling Goff is very close to the gravitational constant G. Then the matter perturbation evolves as Sm oc t^^^ 
during the matter dominance. Meanwhile in the regime m k'^/{a^R) > 1 one has Gofr ~ 4G/(3/,_r), so that the 
evolution of 5m during the matter era is given by 5m o c t(v^-i)/6 [4i|^|4ig. The transition from the former regime 
to the latter regime occurs at the critical redshift |44Cll | 



k 



aoHo 



2n(2n+l) (2(1 - f^D) 



(0)\\2n+l 



,2n 



(f^W)2(„+l) 



l/(6n+4) 



(139) 



where "0" represents the values today. The time at the transition has a scale-dependence t^. oc fc-3/(6n+4)^ which 
means that the transition occurs earlier for larger k. The matter power spectrum Ps^ = |(5n 
acceleration (at time ^a) shows a difference compared to the case of the ACDM model [4161 



at the onset of cosmic 



ACDM 



V33~ 



V33-5 
CX k 8"+* 



(140) 
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The ratio of the two power spectra today, i.e . Pa „ (tn)/-P/i„^'^'^^(^n) is in general different from Eq. (|140[) . but the 
difference is small for n of the order of unity 417 1 . 

The modified evolution o f pe r turbation s for the redshift z < Zk gives rise to the integrated Sachs- Wolfe (ISW) 
effect in CMB anisotropies |42Cll . I44ll4443 |. but this is limited to very large scales (low multipoles). Since the CMB 
spectrum on the scales relevant to the large-scale structure (k > 0.01 /iMpc^^) is hardly affected by this modification, 
there is a difference between the spectral indices of the CMB spectrum and the galaxy power spectrum; Arig = 
(-\/33 — 5)/(6n-|-4). Observationally we do not find any strong signature for the difference of slopes of the two spectra. 
If we take the mild bound An^ < 0.05, we obtain the constraint n > 2. 

The growth index 7 defined in Eq. ([78]) can be as close as 0.4 today in viable f{R) models given in Eqs. (I123p - (|125p 



42l|. Depending on the epoch at which the perturbations cross the transition redshift z^,, the spatial dispersion 



of the growth index 70 = 7(2 — 0) appears in the range of values 0. 40 < 7 < 0.55. There are also regions in 
parameter space for which 70 converges to values around 0.40 < 7 < 0.43 j440l |. These unusual dispersed or converged 
spectra will be useful to distinguish between /(i?) gravity models and the ACDM model in future high-precision 
observations. Since the modified evolution of matter perturbations directly affects the shear power spectrum in weak 
lensing, this is another important test for probing f{R) gravity observationally '444-' 448'| . We also note that the 
nonlinear evolution of matter perturbations in f{R) gravity (corresponding to the scales k > 0.1 /iMpc~^) has been 
studied in Refs. [449l - l457l| . 



3. Local gravity constraints 



Let us discuss local gravity constraints on /(i?) dark energy models. In the region of high density where gravitational 
experiments are carried out, the linear expansion of R in terms of the cosmological value R'^^ and the perturbation 
SR is no lon ger valid because of the violation of the condition SR <C R^^\ In such a nonlinear regime, the chameleon 
mechanism l350l l35l| can be at work to suppress the effective coupling between dark energy and non-relativistic 



matter. In order to study how the chameleon mechanism works in /(i?) grav ity, we transform the action (I116P to the 
Einstein frame action under the conformal transformation g^^ = Fdnf ■ |458| 



Se 



1 

2^ 



1 



(141) 



where K(j) = a/3/2 InF, V{4>) = {RF — /)/(2k^F^), and a tilde represents quantities in the Einstein frame. 

The action (I14ip is the same as d^H]) with the correspondence that in the Jordan frame is equivalent to g'jjil in 
the action (|90)) . Since the quantity F is given by F = e^^'^'^'^ with Q = — 1/\/6 in metric f{R) gravity, the field (p 
is coupled to non-relativistic matter (including baryons and dark matter) with a universal coupling Q ~ —1/^/6. Let 
us consider the models ()123|) and (jl24|) . which behave as Eq. (jl26|) in the region of high density (i? ^ Rc). For the 
functional form p26p the effective potential defined in Eq. ([M)) is 



K=ff(0) 



4k0/\/6 



1 - (2n-|- 1) 



\VGnfiJ 



2n/(2n+l) 



-I- pe 



where 



(142) 



(143) 



Inside and outside a spherically symmetric body, the effective potential (|142p has the following minima given, respec- 
tively, by 



Rc 



2n+l 



K(pB 



Rc 



k^Pb 



2n+l 



(144) 



One has \4>b \ ^ \4>a\ provided that pA ^ Pb- 
The bound (jlOOp translates into 



np. 



^) <1.5xl0-i^ 

PB. 



(145) 



where xi = Ri/R^ and i?i is the Ricci scalar at the de Sitter point. Let us consider the case in which the Lagr angian 
density is given by (|126p for R > Ri. In the original models of Hu and Sawicki [415| and Starobinsky [416l] there 
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are some modification to the estimation of i?i , but this change is not significant when we place constraints on model 
parameters. The de Sitter point for the model (|126p corresponds to fi — 2:^"^^/[2(xf" — n — 1)]. Substituting this 
relation into Eq. (|145p . we find 



2(a;f -n-1) \pB 



n \ 2n+l 

Hi 



<1.5xl0"i^ (146) 



The stability of the de Sitter point requires that m{Ri) < 1, which translates into the condition x^" > 2n^ + 3n + 1. 
Then the term n/[2(a;^" — n — 1)] is smaller than 0.25 for n > 0. Using the approximation that Ri and pB are of 
the orders of the present cosmological density 10~^^ g/cm"^ and the baryo nic / dark matter density 10^^** g/cm^ in our 
galaxy, respectively, we obtain the following constraint from (|146p : [459l | 

n>0.9. (147) 

Thus n does not need to be much larger than unity. Under the condition (|147p . the deviation from the ACDM becomes 
important as R decreases to the order of Re- 

From (|143p we find that there is a curvature singularity with i? -> oo (and — >■ oo) at (/) = for the models 
(|123l) and (I124p . At this singularity the field potential is finite, while its derivative goes to infinity. This singularity 
can be accessible as we go back to the past 460;] , unless the oscillating mode of the scalar- field degree of freedom is 
suppressed. This amounts to the fine-tuning of initi al co nditions for the field perturbation ^416i] . This past singularity 
can be cured by taking into account t he R^ term |46l| . The model of the type f{R) = R — aRc ln(l -I- R/Rc) was 
also proposed to address this problem '462|, but it satisfies neither local gravity constraints [46j| nor observational 
constraints of large-scale structure [464] . Frolov [460] anticipated that the curvat ure s i ngul arity may be accessed in 
a strong gravitational background such as neutron stars. Kobayashi and Maeda j465l 14661 ] showed the difficulty of 
obtaining static spherically symmetric solutions because of the presence of the singularity. On the other hand, the 
choice of accur ate bound ary conditions confirms the existence of static solutions in a strong gravitational background 
with $c < 0.3 [467l - l470l| . 



4- Palatini f{R) gravity 

In the so-called Palatin i formalis m of f{R) gravity, the connections F^^^ are treated as independent variables when 



we vary the action (|116p [47114477 ] . Variation of the action (I116P with respect to g^^, gives 



F{R)R^,{r) ~ -f{R)g^, = , 



(148) 



where F{R) = df/dR, i?^,y(F) is the Ricci tensor corresponding to the connections F^^,, and T^^ is the energy- 
momentum tensor of matter^. Rf^^{T) is in general different from the Ricci tensor calculated in terms of metric 
connections R^^{g). Taking the trace of Eq. (|148p . we find 



F{R)R-2f{R) = K^T, 



(149) 



where T — g^T^^. The trace T directly dete rmin es the Ricci scalar R{T), which is related with the Ricci scalar 
R{9) ^ 9^'^ Rfj-i^ig) the metric formalism via |482j ] 



i?(T) = R{g) 



2(/'(i?(T)))^ 



,(V^/'(i?(T)))(V^/'(i?(r))) 



f'(i?(r)) 



□/'(i?(T)), 



(150) 



where a prime represents a derivative in terms of R{T). Variation of the action (|116p with respect to the connection 
leads to the following equation 



Rf^uig) - -^g,,uR{9) = 



K^T^,, FRiT) - / 



F 
3 

'2F2 



2F 9t.u+p{y^yuF-g^^nF) 



d^Fd,F^-g^,{VFf 



(151) 



* There is another way for the variation of the action, known as the metric-affi ne formah sm [478| |. in which the matter Lagrangian Cm 
depends not only on the metric but also on the connection F^^, . See Refs. [479144831 for the detail of such an approach. 
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Unlike the trace equation (jll8l) in the metric formahsm, the kinetic term OF is not present in the corresponding 
equation (|149p in the Palatini formahsm. Since the time-derivatives of the scalar-field degree of freedom do not appear 
in Palatini f{R) gravity, cosmological solutions are not plagued by the dominance of the oscillating mode in the past. 
In fact, the sequence of radiation, matter, and de Sitter epochs can be realized even for the model f{R) = R — a/R" 
{n > 0) ' 483.-489] . The combined data analysis of SN la, BAO, and the CMB shift parameter places the bound 
n e [-0.23,0.42] on the model f{R) = R-a/R^ {n > -1) [489} . 

Although the background cosmology is well behaved in Palatini f{R) gravity, the evolution of non-relativistic 
matter perturbations exhibits a distinguished feature relative to that in the ACDM model [435, 490 492). Under the 
quasi-static approximation on sub-horizon scales, the equation of matter perturbations is given by |435l | 



2HSr, 



2F 



1 — m 



0. (152) 



Although the matter perturbation evolves as S„i oc t^/^ in the regime |m| k^/{a'^R) ^ 1, the evolution of 6m hi 
the regime \m\ k'^ / {a^ R) ^ 1 is completely different from that in GR. After the perturbations enter the regime 
\m\ k'^ / {a^ R) > 1, they exhibit violent growth or damped oscillations depending on the signs of m |435j |. The 
f{R) models are consistent with observations of large-scale structure if the perturbations do not enter the regime 
|m| fc^/(a^i?) > 1 by today. This translates into the condition 

\m{z = 0)1 < (aoHo/k)'^. (153) 

If we take the maximum wavenumber fc sa 0.2ft,Mpc^"'^ (i.e., k « 600aoi?o)i Eq. (|153p gives the bound \m{z = 0)| < 3x 
10~^. Hence the f{R) models in the Palatini formalism are hardly distinguishable from the ACDM model. 

There are also a number of problems in Palatini f{R) dark energy models associated with the non-dynamical nature 
of the scalar-field d egree of freedom. T he dark energy model f{R) = R — A*^/-R is in confiict with the Standard Model 
of particle physics [475l - l477l 1493^497 1 because of large non-perturbative corrections to the matter Lagrangian. If we 



consi der the models f{R) = R — n"^^"^^^ / R"' , the only way to make such corrections small is to choose n very close 
to |433j |. Hence the deviation from the ACDM model needs to be very small. It was also shown that, for f{R) dark 
energy models, a divergent behavior arises for the Ricci scalar at the surface of a static spherically symmetric star 
with a po lytropic equation of state P = cp^ (3/2 < F < 2), where P is the pressure and po is the rest- mass density 
[497ll498j . These results show that Palatini f{R) dark energy models are difficult to be compatible with observational 
and experimental constraints, although this may not be the case for f{R) models close to the Planck scale [499l - l501 |. 



B. Gauss-Bonnet dark energy models 



It is possible to extend f{R) gravity to more general theories in which the Lagrangian density / is an arbitrary 
function of R, P = R gi^R^" , and Q = RfiuapR^"'^^ , where and R^^ap are Ricci tensor and Riemann tensor 
respectively |502| 5031 . The appearance of spurious spin-2 ghosts can be avoided by taking a Gauss-Bonnet (GB) 
combination f 504| - 507t 



g = R^- AR^.R^" + R^^.c.pR'''"'^ ■ 
A simple dark energy model motivated from low-energy effective string theory [242| is given by [501 



S 



d^x^ 



1 

2^2 



:R 



1 



Sm , 



(154) 



(155) 



where V{(j)) and /(</>) are functions of a scalar field 0, and S]\f is a matt er action. The coupling f{4))Q allows the 
presence of a de Sitter solution even for a runaway field potential V{(j)) '509'-'512l. For the exponential potential 
V{4)) = Voe~'^^'^ and the coupling /(0) = {fo/ p)e'^'^'^ , it was shown in Refs. [509. ,510.] that a scaling matter era can 
be followed by the late-time de Sitter solution ioi p, > \. 

Koivisto and Mota [509] found that the parameter A in this model is constrained to be 3.5 < A < 4.5 (95 % 
confidence level) from the observational data of SN la and WMAP 3-ye ar. The parameter A is constrained to be 
3.5 < A < 4.5 at the 95% confidence level. In the second paper |51ll j. they showed that the model is strongly 
disfavored from the combined data analysis including the constraints coming from BBN, LSS, BAO, and solar system 
data. It was shown in Refs. [510l |513| | that, when the G B term dominates the dynamics, tensor perturbations are 
subject to negative instabilities. Amendola et al. 514| studied local gravity constraints on the above model and 



showed that the energy contribution coming from the GB term needs to be strongly suppressed for the consistency 
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with solar-system experiments. The above results imply that the GB term with the scalar-field coupling f{4')Q can 
hardly be the source for dark energy. 

The d ark energ y models based on the Lagrangian density C — R/{2k^) -{- f{G), have been studied by a number of 
auth ors [5T5l -[52l. In the presence of a perfect fluid with an energy density pM, the Friedmann equation is given by 

[Ullsiel 

3i?2 = gf,g -f- 24H^lggg + PM . (156) 

The Hubble parameter H = Hi dX the de Sitter point satisfies 3i/f = Gif,g{Gi) — f{Gi), where Qi = 24:Hf. The 
stability of the de Sitter point requires the condition < Hff gg{Hi) < 1/384 520]. In order to avoid the instability 
of solutions during radiation and matter eras, we also need the condition f^gg > 0. In Ref. [52d | the authors presented 
a number of f{Q) models that are cosmologically viable at the background level. One of such viable models is given 

by 

f{g) = A-^arctan ( -J ) - aX^/g~, , (157) 



where a, A and are constants. This model can satisfy solar system constraints for a wide range of parameter space 

If we consider cosmological perturbations, however, there is a Ultra- Violet ( UV) instab ility in /{G) models associated 
with a negative propagation speed squared of a scalar- field degree of freedom (516l l524j . This growth of perturbations 
gets stronger on smaller scales, which is difficult to be compatible with the observed galaxy spectrum unless the 
deviation from GR is very small. Thus f{G) dark energy models are effectively ruled out as an alternative to the 
ACDM model. 



Scalar-tensor theories 



There is another class of modified gravity called scalar-tensor theories in w hich the Ricci scalar R is coupled to a 
scalar field ip. One of the simplest examples is Brans-Dicke (BD) theory [29 1| with the action 



S = J d^Xy^ 



Ml 



(158) 



where wbd is the BD parameter, U{ip) is the field pote ntial, and Sm is a matter action that depends on the metric 
gf^i, and matter fields 4'™,. The original BD theory (29ll | does not have the field potential U{(p). 
The general action for scalar-tensor theories can be written as 



S 



d'*x^ 



i/(^,i?)-ic.(^)(V^)2 



+ >5'A/(g^zy, *j\/) , 



(159) 



where / is a general function of the scalar field (p and the Ricci scalar R, u is a function of ip. We choose the unit 
= 1. We consider theories of the type 

f{^,R)^F{^)R-2U{^). 
Under the conformal transformation — Fgf^^^ the action in the Einstein frame is given by [458j | 



(160) 



Se — 



i7?-i(V0)2~y(0) 



where V = U / F"^ . We have introduced a new scalar field 4> in order to make the field kinetic term canonical: 



F 



We define the coupling between dark energy and non-relativistic matter, as 



Q 



Ea 

2F 



F 



F 



3 fF\^ 
F 



-1-1/2 



UJ 

F 



(161) 



(162) 



(163) 
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In f{R) gravity we have lu = and hence F = exp( 1/2/30) from Eq. (|162[) . Then the couphng is given by Q - 
from Eq. ()163|) . If Q is constant as in f{R) gravity, the following relations hold from Eqs. (|162p and (|163|) : 



-■2Q4> 



dip 



In this case the action (|159p in the Jordan frame reads j292 



S 



-F{4>)R - -(1 - 6Q2)f(</,)(V0)2 - t/(0) 



+ SM{g^iu, *m) 



(164) 
(165) 



In the limit that Q the action (|165p reduces to the one for a minimally coupled scalar field with the potential 



-2Q0, 



C/(0). The transformation of the Jordan frame action (|165p under the conformal transformation g^j^^, = e "'^'^g^i, gives 
rise to the Einstein frame action (I16ip with a constant coupling Q. 

One can compare (|165p with the action (|158p in BD theory. S etting = F = e^^*^"^, one finds that two actions are 



equivalent if the parameter wbd is related to Q via the relation [292|, |351 1 



3 + 2wbd = 



2g2 



(166) 



Using this relation, we find that the General Relativistic limit (wbd — ^ oo) corresponds to the vanishing coupling 
(Q — > 0). Since Q — — 1/-\/6 in /(i?) gravity, this corresponds to the Brans-Dicke parameter wbd = j40ll .i525. 526.]. 
One can show tha t the field equation (|151l) in Palatini f{R) gravity is equivalent to the one derived in BD theory 
with wbd = —3/2 |433l 1482^ . Hence Palatini /(i?) gravity corresponds to the infinite coupling (Q^ — >■ oo). 

There are also other scalar-tensor theories that give rise to field-dependent couplings Q(0). For a nonminimally 
coupled scalar field with F[ip) = 1 — ^(^^ and w((/5) = 1 in the action (|159p with p60p . the coupling is field-dependent, 
i.e. Q{ip) = ~ C'/'^(l ~ 6^)]^^^- The cosmological dynamics of dark energy models based on such theories have 

been studied by a number of authors |286l . I527l4532l | . 

Let us consider BD theory with the action (|165[) . In the absence of th e potential U{4)) the BD parameter wbd is 
constrained to be wbd > 4.0 x IC* from solar-system experiments j354l |. This bound also applies to the case of a 
nearly massless field with the potential U{(f) in which the Yukawa correction e"*^'' is close to unity (where M is the 
scalar field mass and r is an interaction length). Using the bound wbd > 4.0 x 10** in Eq. (|166p . we find 



\Q\ < 2.5 X 10 



-3 



(167) 



In this case the cosmological evolution for such theories is hardly distinguishable from the Q = case. Even for 
scalar-tensor theories with such small couplings, i t was sho wn that the phantom equation state of dark energy can be 
realized without the appearance of a ghost state [533l4536l |. 

In the presence of the field potential it is possible for large coupling models {\Q\ = 0(1)) to satisfy local gravity 
constraints under the chameleon mechanism, provided that the mass M of the field 4> is sufficiently large in the 
region of high density. In metric f{R) gravity [Q = — 1/\/6) the field potential ?/((/)) in Eq. (|165|) corresponds to 
U = [FR - /)/2 with (/) = 1/3/2 InF. The viable f{R) dark energy models p23)) and (IT241) have the asymptotic 
form ()126p . in which case the field potential is given by 



2n-l- 1 



(168) 



(169) 



(2n^)2»/(2" 

For BD theories with the constant coupling Q, one can generalize the potential (|168p to the form 

U{(t)) = C/o [1 - C(l - e-2W)P] (C/o > 0, C > 0, < p < 1) . 

As (/) — > 0, the potential (|169p approaches the finite value Uq with a divergence of the field mass squared M"^ = 
U^(f,(j, — )■ 00. This model has a curvature singularity at = as in the case of the f{R) models ()123p and ()124p . The 
mass M decreases as the field evolves away from cf) — 0. The late-time cosmic acceleration can be realized by the 
potential (|169p provided that Uo is of the order of Hq. 

Since the action (I16ip in the Einstein frame is equivalent to the action (|90p . the chameleon mechanism can be 
at work even for BD theories with large couplings {\Q\ = 0(1)). Considering a spherically symmetric body with 
homogenous densities pA and ps inside and outside bodies respectively, the effective potential V^ff — U(0) -I- e'^'^p in 
the Einstein frame (where U(0) = U{(j))/F'^) has two minima characterized by 



J_ f 2UopC 
2Q V PA 



2Q V PB 



(170) 
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Using the experimental bound (jlOOp coming from the violation of e quiva lence principle together with the condition 
for realizing the cosmic acceleration today, we obtain the constraint [292| 

When IQI — 10~^ and \Q\ — 10^^ we have p > 0.68 and p > 0.66, respectively. In f{R) gravity the above bound 
corresponds to p > 0.65, which translates into n > 0.9 for the model p26p . 

The evolution of cosmological perturbations in scalar-tensor theories has been discussed in Refs. j292l . l423l 1434 IsSTI 

l538j |. Under the quasi-static approximation on s ub-horizon scales, the matter perturbation Sm for the theory (|165p 
obeys the following equation of motion |292l l538j 



S„, + 2HS„, - A^^G^sp„Ar^ - , (172) 
where the effective (cosmological) gravitational coupling is 

G (fcVa^)(l + 2Q^)f + 

Here Af^ = t/^^ is the field mass squared. In the "General Relativistic" regime characterized by K'P /F ^ k'^/a?, 
one has Goff ~ G/F and 6m. oc t^f^. In the "scalar-tensor" regime characterized by M"^ / F <^ k'^/a?, it follows that 
Goff ~ (1 -I- 2Q^)GI F and 5m oc f(\/25+48Q^-i)/6^ ^^i^ transition from the former regime to the latter regime occurs 
during the matter era, this gives rise to a difference between the spec tral indices of the matter power spectrum and 
of the CMB spectrum on the scales 0.01 /iMpc"^ <k< 0.2 /iMpc"^ [29^ : 

^^^^il^M,J^I^^^_ (174) 

4 — p 

Under the criterion Arts < 0.05, we obtain the bounds p > 0.957 for Q = 1 and p > 0.855 for Q = 0.5. As long as p 
is close to 1, the model can be consistent with both cosmological and local gravity constraints. 

For the perturbed metric ds^ = —(1 -I- 2'i')dt^ + a'^{t){l — 2^)Siidx^dx^ , the gravitational potentials obey the 
following equations under a quasi-static approximation on sub-horizon scales j292j | 

k^ 4^G ikya^)il + 2Q^)F + AP 

^ ? /L2/„2^^ , n.r2 Z'™'^™ ' (175) 



a 



F {k2/a^)F + AP 



k\ 4^G {kya^)il~2Q^)F + AP 

where we have recovered the gravitational constant G. The results (|175|) and (|176|) include those in f{R) gravity 
by setting Q = — 1/V6. In the regime M^/F <^ k'^/a? the evolution of ^ and $ is subject to change compared to 
that in the GR regi me ch aracterized by M^/F ^ fc^/a^. In general the difference from GR may be quantified by the 
parameters q and C j444| : 

fc2 $ - 

-$ = -47rGgp™<S™ , — — = C ■ (177) 

In the regime M^/F < k'^/a^ of scalar-tensor theory ([TC5]) it follows that 9 ~ (1 - 2Q'^)lF and C -4Q2/(i _ 2Q'^). 

In order to confront dark energy models with the observations of weak lensing, it may be convenient to introduce 
the following quantity [444) 

I] = g(l-C/2). (178) 
From the definition (|177p we find that the weak lensing potential ^ = $ + 'J can be expressed as 

^ = -^TlG'^PmSm^ ■ (179) 

In scalar-tensor theory (|165p one has E = 1/F. The effect of modified gravity theories manifests itself in weak lensing 
observations in at least two ways. One is the multiplication of the term S on the r.h.s. of Eq. (|179l) . Another is the 
modification of the evolution of J,„. The latter depends on two parameters q and ^, or equivalently, E and C- Thus 
two parameters ( E, C) will be useful to detect signatures of modified gravity theories from future surveys of weak 
lensing. See Refs. [539l - l55lj for related works about testing gravitational theories in weak lensing observations. 



34 



D. DGP model 



In the so-caled Dvali, Gabadadze, and Porrati (DGP) [58| braneworld it is possible to realize a "self accelerating 
Universe" even in the absence of dark energy. In braneworlds standard model particles are confined on a 3-dimensional 
(3D) brane embedded in the 5-dimensional bulk space-time with large extra dimensions. In the DGP braneworld model 
[53] the 3-brane is embedded in a Minkowski bulk space-time with infinitely large extra dimensions. Newton gravity 
can be recovered by adding a 4D Einstein-Hilbert action sourced by the brane curvature to the 5D action. Such a 
4D term may be induced by quantum corrections coming from the bulk gravity and its coupling with matter on the 
brane. In the DGP model the standard 4D gravity is recovered for small distances, whereas the effect from the 5D 
gravity manifests itself for large distances. The late-time cosmic acceleration can be realized without introducing a 
dark energy component [552i, ,553 ] (see also Ref. (554] for a generalized version of the DGP model) . 

The action for the DGP model is given by 

I d?X^gR+-^ f d^XV^R- f d^X^CM, (180) 



where qab is the metric in the 5D bulk and g^^i, — d^X^di,X^gAB is the induced metric on the brane with X^{x'^) 
being the coordinates of an event on the brane labelled by x'^. The 5D and 4D (reduced) gravitational constants, 
K^gj and K^^^, are related with the 5D and 4D Planck masses, Af^s) and Af(4), via k^^^ = and k^^^ = l/M'^^y 

The first and second terms in Eq. (|180|) correspond to Einstein-Hilbert actions in the 5D bulk and on the brane, 

rane\ 



respectively. The matter action consists of a brane-localized matter whose action is given by J d^x^— g (it + £ 
where a is the 3-brane tension and C^f'^'' is the Lagrangian density on the brane. Since the tension is not related to 
the Ricci scalar R, it can be adjusted to be zero. 

The Einstein equation in the 5D bulk is given by G\g = 0, where G\g is the 5D Einstein tensor. Imposing the 
Israel junction conditions on the brane with a Z2 symmetry, we obtain the 4D Einstein equation ^b. 



Gf,^ - — {Kf,^ - Xg,,^) = K(4)?)„. , (181) 

where Kf^i, is the extrinsic curvature on the brane and Tf^u is the energy-momentum tensor of localized matter. The 
cross-over scale rc is defined by Tc = k^^^/ {2k^^.^). The Friedmann equation on the fiat FLRW brane takes a simple 



form 1552 ^ 

,2 



e 



K' 



H'~-H^ -^PM , (182) 

Tc 3 

where e = ±1, and pM is the energy density of matter on the brane (with pressure Pm) satisfying the continuity 
equation 

Pm + 3H{pm + Pm) = 0. (183) 

If Tc is much larger than the Hubble radius H^^, the first term in Eq. (I182p dominates over the second one. In this case 
the standard Priedmann equation, — k^^.^pm/3, is recovered. Meanwhile, in the regime rc < H^^, the presence of 
the second term in Eq. (|182l) leads to a modification to the standard Friedmann equation. In the Universe dominated 
by non-relativistic matter (pM oc a~^), the Universe approaches a de Sitter solution for e ~ +1: H — >■ H^s = l/j^c- 
Hence it is possible to realize the present cosmic acceleration provided that rc is of the order of the present Hubble 
radius Hq^. 

Although the DGP braneworld is an attractive model allowing a self acce leration, the joint constraints from SNLS, 
BAO, and CMB data shows that this model is disfavored observationally 5561456 1| . There is a modified version of 



the DGP model characterized by the Friedmann equation H'^ — H" /r1 " = K^^^pml'^, where is the energy density 



)(0) 



of non-relativistic matter 562| . In Fig. IHlwe show la and 2a contours in the (Om , a) plane constrained from the 
joint data analysis of SN la, BAO, CMB, gamma ray bursts, and the linear growth factor of matter perturbations 
[56ll |. The parameter a is constrained to be a = 0.254 ± 0.153 (68 % confidence level) and hence the flat DGP model 
{a = 1) is incompatible with observations. 

The evolution of density perturbations in the DGP model has been studied in Refs. (56314571 ]. Under the quasi- 



static approximation on sub-horizon scales, the linear matter perturbation 5m (with a homogenous density pm) obeys 
the following equation 

^™ + 2i/<5„, - 47rGefrp™<5™ = , Ges = [1 + 1/(3/?)] G , (184) 
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Figure 6: Combined observational constraints on the modified DGP model characterized by the Friedmann equation — 
i^^/r^"" = K(4)Pm/3. The two curves show la and 2(t contours in the ,ct) plane. The original DGP model (q = 1) is 
incompatible with observations. From Ref. 561]. 



where /3 = 1 — 2Hrc[i + H/{3H'^)]. In the deep matter era {Hrc ^ 1), /3 is largely negative ^ 1). In this 
regime the matter perturbation evolves as Sm oc t^^^ as in GR. Around the late-time de Sitter solution one has 
/3 ~ 1 — 2Hrc — —1 and 1 + 1/(3/3) ~ 2/3, so that the growt h rat e of 5m gets smaller than that in the ACDM model. 
The growth index 7 defined in Eq. ([75)1 is given by 7 « 0.68 323 1. 

In the massless regime {M'^/F ^ k'^/a?) in BD theory the effective gravitational coupling (|173l) is given by 
Goft = (1 + 2Q'^)G/F = (4 + 2wbd)/(3 + 2lobt>)G/F, where we used the relation (ITBS| . Comparing this with the 
effective coupling (|184l) . we find that the DGP model is related to BD theory via wbd = (3/2)(/? — 1) with F ^ 1. 
Since /3 < for the self- accelerating DGP solution, this implies that wbd < — 3/2 and hence the DGP model contains 
a ghost mode. It is however possible to construct a generalized DGP model free from the ghost problem by embedding 
our visible 3-brane with a 4-brane in a flat 6D bulk 168]. 

In the DGP model, a brane bending mode in the bulk corresponds to a scalar-field degree of freedom. In general 
such a field can mediate a long range fifth force incompatible with local gravity experiments, but the presence of a 
self-interaction of </> allows the so-called Vainshtein mechanism |572| to work within a radius = {rgr^Y/^ {vg is the 
Schwarzschild radius of a source) . The DGP m odel can be consistent with local gravity constraints under some range 
of conditions on the energy-momentum tensor [573l - [575t . 

The DGP model stimulated other approaches for constructing ghost-free theories in the presence of nonlinear self- 
interactions of a scalar field 0. It is important to keep the field equations at second order in tim e derivatives to avoid 
that an extra degree of freedom gives rise to a ghost state. In particular Nicolis et al. [577j ] imposed a constant 
gradient-shift symmetry ("Galilean" symmetry), df^<j} — )■ d^(j) + b^, to restrict the equations of motion at second order, 
while keeping a universal gravitational coupling with matter. In the 4-dimensional Minkowski space-time they found 
five terms Ci (z = 1, • • • 5) giving rise to equations of motion satisfying the Galilean symmetry. The first three terms 



are given by £1 — 4>, C2 = V^(/)V'^0, and £3 = □(/(V^^V' 



appears in the DGP model, which allows the possib ility for the consistency with solar system experiments through 
the Vainshtein mechanism. Deffayet et al. |578l l579l ] derived the covariant expression of the terms £i {i 



The term £3 is the nonlinear field derivative that 

rough 
•5) by 



1, 



extending the analysis to the curved space-time. 

Silva and Koyama [5801] considered BD theory (without a field potential) in the presence of a nonlinear derivative 
term ^(0)n0V^(/)V''(/), i.e. 



^((/.)n0v^(/.v^ 



s 



M 



(185) 



where Sm is the matter action. Although the nonlinear derivative term does not satisfy the Galilean i nvari ance in 
the FLRW cosmological background, the field equation of motion remains at second order (see also Ref. j58l| ). This 
term also arises as one of higher-order derivative corrections to low-energy effective string theory. Interestingly, 
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for the function ^(0) = 1/(M^(/)^), there exists a de Sitter solution responsible for dark energy (provided that 
M « Hq). Moreover, because of the presence of the nonlinear interaction, t he pr oblems of the appearance of ghosts 
and instabilities can be avoided for the BD parameter ujbd smaller than —2 [580j |. At early times General Relativity 
can be recovered by the cosmological Vainshtein mechanism. A number of interesting observational signatures such as 
modified growth of matter perturbations as well as a distinguished ISW effect have been studied in Refs. [58Cki582i . i583l] . 



VI. COSMIC ACCELERATION WITHOUT DARK ENERGY 



There are attempts to explain the apparent cosmic acceleration by inhomogeneities in the distribution of matter 
without recourse to a dark energy component (see Ref. [584l | for review). 

One of such approaches is the void model in which the presence of underdense bubbles leads to the faster expansion 
of the Universe compared to the outside. In other words we live in the middle of a huge spherical region and we 
interpret the evolution of this underdense region as an apparent cosmic acceleration. Originally, Tomita [s^ HO] 
introduced a local homogenous void separated from the outside described by a homogenous FLRW space with a 
singular mass shell. The analysis was extended to the models with a continuous transition between the inside and 
outside the void [i^. This can be described by a class of the Lemaitre-Tolman-Bondi (LTB) spherically symmetric 
models. Theoretical and observational aspects of the LTB model have been extensively studied as an alternative to 
dark energy [61-63, 585l-[610t. 

The second approach is based on the backreaction of cosmological perturbations arising from perturbing the ho- 
mogeneous Universe (63 - [66| . Unlike the void model, this tries to explain a real cosmic acceleration by arranging 
inhomogeneities that come from the deviation of the FLRW metric. 

There is another a pproach fo r explaining the apparent cosmic acceleration based on the "Ultra Strong" version of 
equivalence principle |61ll |612| | . In this model the standard geometric description of space-time as a metric manifold 
holds as a small distance approximation and hence General Relativity can be m odified on large scales by a curvature- 
dependent subleading effect. Although the original model proposed in Ref. [613j | do not explain the observational data 
of SN la very well, its modified version can be consistent with the SN la data with a rather low value of the Hubble 
constant, Hq « 50kmsec~^ Mpc^^. 

In the following we shall briefiy review the first two approaches. 



A. Inhomogeneous Lemaitre-Tolman-Bondi model 

In order to discuss a spherical inhomogeneity in local regions we take the LTB metric given by 

ds^ = -dt^ + X^{t, r) dr^ + R^{t, r) dQ^ , (186) 

where the expansion factor along the radial coordinate r is different relative to the surface line element dft^ — 
d9^ + siia^ 9d(jP. Solving the (0,1) component of the Einstein equation Gqi = for the fiuid at rest, it follows that 
X{t,r) is separable as X{t,r) = R\t , r) / ^J\ + I3{r) . Here a prime represents a partial derivative with respect to r 
and (5{r) is a function of r. Then the metric (|186p is given by 

= -dt^ + M^dr2 + R\t, T)dif . (187) 

1 + I3(r) 

The metric p87p recovers the one in the FLRW space-time by the choice R = a{t)r and /3 = —Kr^, where K is 
a cosmic curvature. In other cases the metric (11871) describes a spherical inhomogeneity centered on the origin. We 
define the transverse Hubble function H± and the radial Hubble function as H± = R' /R' and H± — R/R. The 
Einstein equations in the presence of non-relativistic matter (energy density pm) give (63j 

i/i + 2i/||i/^ - 4 - 1^ = 8^Gp™ , (188) 



i?2 RR 



6- + 2Hi-2^-2H^^H^ + -^ = -87rGp,„ . (189) 



Eliminating the term pm from Eqs. (|188p and (|189p . we obtain the relation 2RR + R^ — /3{r). Integrating this 
equation, it follows that 



(190) 



37 



where a(r) is an arbitrary function of r. From this we can introduce the today's effective density parameters of matter 
and the spatial curvature, respectively, as 

- ' ir) = l- nl^^ ir) = ^ . (191) 

We define the time t = at the decoupling epoch (the redshift z ~ 1090) with R{r,t = 0) = 0. Introducing the 
conformal time 77 as d?? = {s/^/R) dt, we obtain the following parametric solutions of Eqs. (|188p and (|189p for /3 > 

"^'^^ -(sinh77-ry) = (sinhTy-r/). (193) 



The structure of the void with an under d ensit y can be accommodated by choosing fim^ (r) and h = 
i/j^o/(100kmsec^^ Mpc^^) in the following form j598l | (see also Ref. for another choice) 

4?(r) = r!out + (r!in-f^out)/(r,ro,A), (194) 
Hr) = /iout + (/iin-/iout)/(r,ro,A), (195) 

where the function /(r, ro, A) = [1 — tanh((r' — ro)/2A)]/[l + tanh(ro/2A)] describes a transition of a shell of radius 

To and thickness A ("in" and "out" represent quantities inside and outside the void, respectively).^ 

The trajectory of photons arriving at r = today is characterized by a path t — i{r) satisfying |590j | 



di _ R'{r,i) 
d^ ~ ~ ^/TTp' 

Then one can show that the redshift z — z{r) of photons obeys the differential equation (62 



(196) 



with z{r = 0) = 0. The luminosity distance dL{z) is related to the diameter distance dA[z) = R{r, i) according to the 
usual duality relation 

dL{z) = {l + zfR{r,t) . (198) 

Now we are ready to confront the inhomogeneous LTB model with the SN la observations. From the requirement 
that the CMB acoustic peak is not spoiled we require that local density parameter fiin is in the range 0.1-0.3, whereas 
f^out = 1 as predicted by inflation. The observed local value of H is around /lin « 0.7, whereas outside the void 
one requires h « 0.5 t o be consistent with i^out — 1- Then the two parameters rp and A are constrained by the SN 
la data. In Ref. ' 598| it was shown that the inhomogeneous LTB model can be consistent with the SN la data for 
ro = 2.3 ± 0.9 Gpc and A/ro > 0.2. 

We note, however, that one can place other constraints on the void model. If we do not live around the center of 
the void, the observed CMB dipole becomes much larger than that allowed by observations. The maximum distance 
Tc to the center is constrained to be smaller than 10-20 Mpc [590 , 605]. Even if we happen to live very close to the 
center of the void, we observe distant off-centered galaxy clusters. Such off-centered clusters should see a large CMB 
dipole in their reference frame. For us this manifests itself observationally as a kinematic Sunyaev-Zeldovich effect. 
Using the observ ational da ta of only 9 clusters, the inhomogeneous LTB model with void sizes greater than 1.5 Gpc 
can be ruled out |599l . |604| . This is already in mild conflict with the constraint derived by the SN la data. It remains 
to see whether the void model can be ruled out or not in future observations. 



B. Backreaction of cosmological perturbations 



Let us finally discuss the possibility of realizing a real cosmic acceleration by the backreaction of inhomogeneities to 
the FLRW space-time. In general averaging the inhomogeneities and then solving the Einstein equations (the standard 
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approach) might not be the same as solving the full inhomogeneous Einstein equations first and then averaging them. 
In other words, the expected value of a nonlinear function of x is not the same as the nonlinear function of the 
expected value of x. The argument is complicated and controversial, so we mention the basic ideas onl y briefly . Th e 
readers who are interested in the detail of this line of research may have a look at the original papers 



Let us decompose the Einstein tensor G^jy and the energy-momentum tensor T^i, into the background (0-th order 

i, as 

component of the Einstein equation ^ gives 



FLRW Universe) and the perturbed parts, as G^i, — G^"J + G^^} and T^i, — T^^J + T^u) , respectively. Then the (00) 



4o^=8^G(t(°)+tW)-G«. (199) 
Identifying the average matter density at this order as {p) = + and averaging over Eq. (|199p . we obtain 

(g|,°)) = 87rG(p) - {G^^) , (200) 

where Gqq'' = 3H^ in the FLRW background. This shows that ^ 87rG(p) in general because of the presence of 

the term (g[,o^). If we first average the metric as (gl^u) = 0, it then follows that Gqq'' = Goo{{gl^u)) = and hence 
= 8nG{p). The above argument can be extended to the second order, i.e. 

{Gi'S)^8nG{p)^{Gg^+Gg^). (201) 

The cosmological evolution depends on the averaging procedure. For example, if we take the average of the function 
f{t,x) as [H 



,,,,^^_/d3./^/(M) (202) 

Jd^X^/J{t,x) 

where 7 corresponds to the determinant of the perturbed metric of spatial constant-time hypersurfaces, then the 

(2) r 

second-order term Gqq contributes to the expansion rate of the Universe (which is typ ically of the order of 10^""). If 
we use other ways of averaging, the amplitude of such a term is subject t o cha nge j616l |. Also the results are affected 
by adding the contributions higher than the second order. In fact Ref. [614| pointed out the danger of arbitrarily 
stopping at some order by showing several examples in which many contributions cancel each other. 

The backreaction scenario is very attractive if it really works, because it is the most economical way of explaining 
the cosmic acceleration without using dark energy. We hope that further progress will be expected in this direction. 



VII. CONCLUSIONS 



We summarize the results presented in this review. 

• The cosmological constant (wde ~ ~1) is favored by a number of observations, but theoretically it is still 
challenging to explain why its energy scale is very small. 

• Quintessence leads to the variation of the field equation of state in the region > —1, but the current 
observations are not sufhcient to distinguish between quintessence potentials. 

• In k-essence it is possible to realize the cosmic acceleration by a field kinetic energy, while avoiding the instability 
problem associated with a phantom field. The k-essence models that aim to solve the coincidence problem 
inevitably leads to the superluminal propagation of the sound speed. 

• In coupled dark energy models there is an upper bound on the strength of the coupling from the observations 
of CMB, large-scale structure and SN la. 

• The generalized Chaplygin gas model allows the unified description of dark energy and dark matter, but it needs 
to be very close to the ACDM model to explain the observed matter power spectrum. There is a class of viable 
unified models of dark energy and dark matter using a purely k-essence field. 

• In f{R) gravity and scalar-tensor theories it is possible to construct viable models that satisfy both cosmolog- 
ical and local gravity constraints. These models leave several interesting observational signatures such as the 
modifications to the matter power spectrum and to the weak lensing spectrum. 
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• The dark energy models based on the Gauss-Bonnet term are in conflict with a number of observations and 
experiments in general and hence they are excluded as an alternative to the ACDM model. 

• The DGP model allows the self-acceleration of the Universe, but it is effectively ruled out from observational 
constraints and the ghost problem. However, some of the extension of works such as Galileon gravity allow the 
possibility for avoiding the ghost problem, while satisfying cosmological and local gravity constraints. 

• The models based on the inhomogeneities in the distribution of matter allow the possibility for explaining the 
apparent accelerated expansion of the Universe. The void model can be consistent with the SN la data, but it 
is still challenging to satisfy all other constraints coming from the GMB and the kinematic Sunyaev-Zeldovich 
effect. 

When the author submitted a review article Q on dark energy to International Journal of Modern Physics D in 
March 2006, we wrote in concluding section that "over 900 papers with the words 'dark energy' in the title have 
appeared on the archives since 1998, and nearly 800 with the words 'cosmological constant' have appeared". Now in 
April 2010, I need to change the sentence to "over 2250 papers with the words 'dark energy' in the title have appeared 
on the archives since 1998, and nearly 1750 with the words 'cosmological constant' have appeared". This means that 
over 4000 papers about dark energy and cosmological constant have been already written, with more than 2300 papers 
over the past 4 years. Many cosmologists, astrophysicists, and particle physicists have extensively worked on this new 
field of research after the first discovery of the cosmic acceleration in 1998. We hope that the future progress of both 
theory and observations will provide some exciting clue to reveal the origin of dark energy. 
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